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SECTION  I 
INTRODUCTION 


As  part  of  Contract  No.  F33615-83-C-3020  the  University  of 
Dayton  has  conducted  a  study  of  technologies  required  for  an 
interferanemography  system.  Interferometry  and  laser  Doppler 
anemometry  were  evaluated  for  their  applicability  to  inter¬ 
feranemography,  as  were  other  technologies. 

The  specified  requirements  for  an  interf eranemographic 
system  are  to  provide  real-time/on-line  simultaneous  measurements 
of  density  and  velocity  in  complex  flow  fields. 

Interferometry  is  a  generally  accepted  method  for  measure¬ 
ment  of  air  density.  Within  the  past  decade  holographic  inter¬ 
ferometric  techniques  have  been  used  successfully  to  obtain 
density  data  in  a  variety  of  wind  tunnel  studies.  However,  the 
general  application  and  wide  use  of  interferometry  is  hindered  in 
wind  tunnel  studies  by  the  following:  (a)  single  line-of-sight 
views  of  the  flow  field  which  limit  the  application  to  2-D  or  3-D 
axisymmetrlc  fields,  (b)  insufficient  measurements  for  optically 
weak  fields  (i.e.,  hypersonic  flows  or  low-density  flows),  (c) 
manual  data  reduction  methods  now  widely  used.  The  application 
of  traditional  Interferometry  to  flow-field  measurements  is  pre¬ 
sented  in  Appendix  A,  which  is  the  Interim  Technical  Report 
UDR-TR-84-48  .  In  view  of  the  limitations  of  standard  inter¬ 
ferometric  techniques,  the  Lau  interferometer  was  adapted  for  use 
in  flow-field  measurements.  This  has  led  to  some  novel  opera¬ 
tional  capability  for  interferometers,  including  measurement  in  a 
restricted  volume. 

Laser  Doppler  techniques  have  been  used  to  measure  veloci¬ 
ties  in  subsonic  and  supersonic  flow  fields.  Laser  Doppler 
Anemometry  (LOA)  is  a  process  whereby  changes  in  the  frequency 
(Doppler  shifting)  of  laser  radiation  scattered  from  particles 
moving  through  an  optical  control  volume  are  monitored  to  compute 
the  velocity  of  the  moving  particles. 
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Appendix  A,  the  Interim  Technical  Report,  presents  a  devel¬ 
opment  of  the  requirements  and  techniques  for  Interferometric 
measurements  of  a  general  three-dimensional  flow  field.  As  a 
part  of  the  study  using  interferometry,  it  was  necessary  to 
characterize  the  interferometer  sensitivities  for  the  different 
wind  tunnels  at  Wright-Patterson  Air  Force  Base,  and  the  results 
of  this  comparison  are  presented  in  Section  2,  Appendix  A. 

Section  3,  Appendix  A,  presents  a  development  of  the  LDA 
requirements  and  techniques  for  characterization  of  local  veloci¬ 
ties  in  a  complex  flow  field.  Alternate  techniques  for  measure¬ 
ment  of  the  local  velocities  are  compared  to  LDA  to  see  if  they 
offer  any  advantage  over  LDA  in  work  by  the  Flight  Dynamics 
Laboratory. 

Section  4,  Appendix  A,  summarizes  the  potential  of  existing 
technology  to  obtain  real-time/on-llne  Interf eranemographic 
measurements  in  complex  flow  fields. 

The  theory  of  the  Lau  interferometer,  as  it  was  used  for 
flow-field  measurement,  is  developed  in  Section  II  of  this  reprt. 
The  system  performance  was  demonstrated  through  visualization 
measurements  of  the  varying  density  fields  of  a  freon  jet  exiting 
into  air  and  of  a  thermal  field  around  a  heated  cylinder.  These 
are  presented  in  Section  III  of  this  report. 


SECTION  II 

INTERFEROMETRIC  MEASUREMENT  OF  LOCALIZED  FLOW-FIELD  DENSITIES 


The  goal  of  this  program  was  to  provide  an  interferometric 
technique  capable  of  making  localized  measurements  of  density  in 
a  flow  field.  This  is  to  be  combined  with  a  velocity  probe.  The 
velocity  probe  analysis  was  submitted  as  an  Interim  Technical 
Report,  and  is  attached  as  Appendix  A.  Included  in  that  interim 
report  was  a  discussion  of  the  status  of  interferometry  and  its 
application  to  density  measurements.  The  basic  problem  is  that 
interferometric  techniques  integrate  the  phase  variation  along 
the  optical  path.  In  this  report  we  will  discuss  our  attempts  to 
find  a  localized  measurement  technique.  The  study  has  led  to  the 
adaptation  of  the  Lau  Interferometer.  The  optical  configuration, 
its  method  of  operation,  and  its  potential  for  localized  measure¬ 
ments  will  be  described  in  the  following  paragraphs. 

1.  INTERFEROMETRY  PRINCIPLES 

As  a  ray  of  light  traverses  a  medium,  the  phase  of  the 
electromagnetic  radiation  at  any  location  is  dependent  upon  the 
refractive  index  of  the  medium  Integrated  along  the  optical 
path.^  If  two  rays  travelling  different  optical  paths  are  mixed, 
the  resulting  fringe  shift  is  given  by 


S(x,y)  =  X  ^  (n(x,y,z)  -  dl 

[see  Appendix  A,  Equation  (1)].  This  equation  shows  that  the 
effect  observed  is  the  result  of  the  difference  in  index  of 
refraction  of  the  test  ray  from  the  reference  ray,  integrated 
along  the  entire  optical  path.  To  limit  the  region  of  obser¬ 
vation  it  is  necessary  to  provide  a  common  path  for  the  interro¬ 
gation  ray  and  the  reference  ray  over  the  entire  path  length, 
except  in  the  region  of  interest.  This  will  yield  a  local 
measurement  of  index  change  which  can  be  related  to  the  local 
density  variation  through  the  Dale-Gladstone  coefficient 


by  the  equation 


Ap  =  An/Kjjg 


2 .  OPTICAL  SYSTEM 

The  Initial  concept  was  based  upon  the  use  of  an  optical 
system  using  fast  lenses  and  a  white-light  source.  A  white-light 
illuminated  grating  was  focused  in  the  region  of  interest  and 
then  reimaged  on  a  matching  grating.  The  transmission  of  the 
combined  gratings  was  then  imaged  onto  a  detector  (Figure  1). 
Since  the  source  was  incoherent,  every  point  on  the  first  grating 
would  be  an  effective  source,  Independent  of  every  other  position 
on  the  grating.  By  appropriately  selecting  components  the 
radiation  from  each  source  would  propagate  along  a  common  optical 
path  for  most  of  the  region.  Only  in  the  region  of  interest 
would  the  optical  paths  be  distinctly  separate.  Once  past  this 
region  the  optical  paths  would  again  merge.  As  a  result  of  the 
common  path,  only  phase  variation  in  the  region  of  interest 
results  in  relative  path  length  variation.  In  the  original  con¬ 
cept  th  phase  variation  would  lead  to  a  deflection  of  the 
emerging  ray  and  a  deflection  of  the  image,  resulting  in  a  moire 
effect.  Unfortunately,  a  lens  is  the  ideal  component  for  re- 
imaging  the  deflected  ray  back  onto  the  same  location  where  the 
undeflected  ray  reimages  and,  as  a  result,  no  effect  was  observed 

The  concept  was  considered  to  be  still  valid;  the  question 
was  how  to  convert  a  deflection  into  a  movement  of  the  image. 

The  answer  comes  from  Fourier  optics.  A  change  in  direction  in 
the  front  focal  plane  of  a  lens  is  equivalent  to  a  deflection  due 
to  a  spatial  index  variation  and  it  will  result  in  a  translation 
of  the  image  in  the  back  focal  plane  by  an  amount  proportional  to 
the  angle  (spatial  frequency).  Of  the  many  potential 


realizations  of  such  a  system  we  studied  two  in  detail  and 
implemented  one  in  the  laboratory. 

The  concept  of  the  system  implemented  is  simple  to 
visualize  when  using  a  coherent  source  and  then  straightforward 
to  extend  to  an  incoherent  source.  When  coherent  light  illumi¬ 
nates  a  subject  in  the  front  focal  plane  of  a  lens,  the  spatial 
Fourier  transform  is  imaged  at  the  rear  focal  plane  of  the  lens. 

If  the  subject  is  a  Ronchl  ruling,  the  spatial  Fourier  transform 
is  a  series  of  equally  spaced  dots,  with  a  well  known  intensity 
variation,  along  a  line  through  the  optical  axis  in  a  direction 
normal  to  the  ruling.  The  second  lens  retransforms  the  rays, 
reforming  the  image  (Figure  2).  If  there  is  a  spatial  phase 
variation  at  the  original  Fourier  plane,  the  resulting  image 
would  vary  in  a  manner  related  to  the  phase  variation.  In 
reforming  the  image  onto  a  matched  ruling,  variations  in  the 
reformed  image  would  become  apparent.  To  overcome  the  presence 
of  dead  spots  in  the  Fourier  plane — spatial  locations  where  there 
is  no  light  to  Interrogate  the  local  density — the  system  was 
changed  from  a  coherent  source  to  an  incoherent  white  light 
source.  The  wide  range  of  wavelengths  results  in  a  uniform 
coverage  at  the  front  focal  plane.  The  validity  of  employing  an 
Incoherent  source  is  discussed  in  Section  III. 

The  other  potential  system  considered  is  closely  related  to 
the  first.  The  optical  system  would  provide  two  optical  beams  at 
a  small  angle  to  each  other  but  overlapping  over  most  of  their 
passage  through  the  optical  system.  They  would  distinctly 
separate  only  near  the  focal  plane.  The  reimaged  system  would  be 
reformed  into  a  single  path,  focused  onto  a  detector,  and  mixed  in 
the  detection  process.  This  system  can  be  implemented  using 
Wollaston  prisms  as  shown  in  Figure  3.  The  observed  signal  would 
be  the  variation  in  intensity.  In  this  format  the  signal  is  a  coarse 
measurement  tool.  Several  optical  systems  with  much  higher 
resolution,  based  on  this  technique,  can  be  realized  in  a  straight¬ 
forward  manner.  For  example,  a  Bragg  cell  used  as  a  beamsplitter 


would  also  introduce  a  frequency  offset  and  produce  a  bias  fre¬ 
quency  making  phase  sensitive  detection  usable.  The  system  might 
be  worth  pursuing  in  a  separate  investigation. 

The  optical  system  selected  in  the  current  study  is  a  Lau 
interferometer.  The  application  of  this  system  to  the  measure¬ 
ment  of  flowfields  is  novel  and  we  will  report  the  results.  The 
system  is  very  robust,  is  easy  to  align,  and  uses  a  bright  white- 
light  source.  All  of  these  properties  are  advantageous  when 
taking  measurements  in  wind  tunnels  where  there  is  often  a  high- 
level  acoustic  environment,  where  it  is  difficult  to  get  good 
optical  access,  and  where  light  collection  efficiency  is  poor. 

3.  THE  LAU  INTERFEROMETER 

It  is  well  known  that  incoherent  light  can  gain  a  measure 
of  coherence  merely  by  propagating  through  space  as  described  by 
the  van  Cittert-Zernlcke  Theorem. 2  we  will  denote  the  state  of 
coherence  by  the  cross-spectral  density  function^  W(x',x”,«) 
which  describes  the  phase  relationship  between  two  points  x'  and 
x"  as  a  function  of  the  frequency  w  or  equivalently  the  wave¬ 
length  X.  The  cross-spectral  density  function  may  be  propagated 
through  a  system  by  a  series  of  Fresnel-like  integrals  in  the 
same  manner  as  the  mutual  coherence  function.*  The  cross- 
spectral  density  function  has  the  added  benefit  that  it  is  a 
valid  description  for  broadband  illumination;  the  more  commonly 
used  mutual  coherence  function  is  valid  only  in  the  quasi-mono- 
chromatic  approximation. 

We  can  use  the  propagation  Integral  to  find  the  relation 
between  the  coherence  of  light  fields  on  opposite  sides  of  a 
lens,  a  situation  common  to  most  optical  systems.  At  the 
distance  d  in  front  of  the  lens  an  object  with  complex  amplitude 
transmittance  G(Xj)  is  illuminated  by  light  with  a  cross-spectral 
density  function  (Xj^ '  ,Xj^" )  (the  u  is  dropped  for  convenience). 
The  light  propagation  to  the  lens  is  represented  by  exp [-irrx,^/XF’j 


We  calculate  the  cross-spectral  density  function  W^Cx'^.x"^) 

In  a  plane  located  F-A  away  from  the  lens.  It  will  be  convenient 
later  to  consider  displacements  about  the  focal  plane,  but  A 
should  not  be  considered  a  small  quantity.  By  repeated  applica¬ 
tion  of  the  double  Fresnel  Integral  we  arrive  at 


=  C  //  e 


-l-^x'2  1-2-  x" 2 

XF  *2  ^XF  *2 


t  II  I  •  II 


e  //w^(Xj^,x^)G(x^)6^(Xj) 


(1) 


I  I  o  *  • 


^Xd  -(Xg-x^)  ]  ^  ^X(F-A)  '•'"3  “^2 

X  e  dx^dx^e 


I  I  O  II  M 


[{X,-X„)  -(Xg-Xj)  ] 

dXgdXg  . 


The  factor  C  contains  the  constants  of  propagation.  It  Is  con¬ 
venient  to  switch  to  a  sum  and  difference  notation.  We  let 

I  If  I  II 

x^  =(x^  +  x^),  Ax^  “  *1  ~  *1’  Then  equation  (1)  becomes: 


W(Xg,AXg)  =  C  JJ  e 


.2n  . 

•‘Xp  *2**2 


AX^  Ax. 

fj  W^(Xj,AXj)G(x^+  ^)G*(X^--^  ) 


^  Xd  [XgAXg+X^Ax^-XjAXg-XgAXj]  1  WF^ 


2ji 


X  e 


X  dx^dAx^dXgdAXg  , 


X(F-A)  i''3“''3"''2“''2  ''2“'‘3  ''3“'‘2- 

(2) 


We  note  that  by  letting  AXg  »  0  we  find  the  Intensity  In  plane  3. 

After  a  great  deal  of  tedious  but  straightforward  manipula¬ 
tion  of  Equation  (2)  we  arrive  at 
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(X,AX  +X,AX  ) 
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where  a  «  1-d/F. 

We  cam  distinguish  several  cases  of  interest.  If  A  »  0 
(observation  in  the  focal  plane  of  the  lens)  then 

if^  w  ^*1  • 

V)^(X^.LX^)  =  Ce  ^Jf  W^(Xj.AXj)6(Xj+  - 


^  (Xi»X3«34x.) 


For  any  state  of  coherence  W^  there  is  a  two-dimensional 
Fourier  transform  relationship  between  and  W^. 

If  Wj(x^,AXj)  =  6(Ax^),  then  the  input  illumination  is 
totally  incoherent.  Equation  (3)  becomes 
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where  G  is  the  Fourier  transform  of  |g|  ,  the  intensity 

transmittance  of  the  object. 

Equation  (5)  is  a  significant  result.  At  all  distances 
beyond  the  lens  the  field  coherence  is  in  the  form  of  a  Fourier 
transform  of  the  input  object.  The  effect  of  changing  A  is  to 


y  v*".-  '■  ’> 


change  the  scaling  of  the  Fourier  transform.  The  coordinate  of 
the  Fourier  transform  is  the  difference  coordinate  Ax.  That  is, 
if  a  point  is  fixed  in  the  observation  plane.  Equation  (5) 
gives  the  coherence  relation  a  distance  Ax  away  from  x^.  There 
is  no  fixed  reference  system  for  evaluating  Equation  (5);  it  must 
be  evaluated  for  every  point  in  the  plane  of  Interest. 

We  will  consider  the  special  case  of  a  grating-like  input 
in  plane  1 .  Let 

00 

G(x)  =  T(x)  *  X  6(x-nb), 

n=— 00 

where  *  represents  convolution.  Equation  (6)  describes  a  func¬ 
tion  T{x)  replicated  about  x=nb,  n=0,  +1,  ±2,  ...  T(x)  can 
be  considered  a  grating  line  shape  and  has  a  width  less  than  b. 

In  the  case  of  a  simple  grating  or  Ronchl  ruling  T(x)  is  either  1 
or  0,  so  T(x)  =  |T(x)|^‘  Equation  (6)  then  becomes 

F-Ao  ^  Ax-  "  Ax 

W3(Xa,AXa)  =  ce  T  I  *^A(F-A)  "  b 

m=— 00 

That  is,  the  coherence  function  is  a  set  of  delta  functions 
spaced  at  Ax^  =  mA(F-A)/b  and  weighted  by  T( Ax^/A (F-A ) ) .  If  a 
point  x^  is  chosen  as  a  reference,  then  a  point  at  x^+A(F-A)/b 
will  be  coherent,  but  a  point  immediately  adjacent  to  it  will 
not . 

The  fact  that  an  Incoherently  illuminated  grating  will  give 
rise  to  a  grating-like  coherence  function  is  the  basis  of  the  Lau 

5 

effect.  A  second  grating  illuminated  by  the  coherence  function 
of  Equation  (7)  will  sample  the  coherence  function.  If  the 
distance  between  gratings  is  properly  chosen,  the  second  grating 
will  only  sample  points  that  are  coherent  with  respect  to  each 


other.  That  light  will  propagate  to  the  far-field  and  form 
interference  fringes.  Thus,  the  Lau  effect  describes  a  coherent 
interference  phenomenon  which  arises  from  incoherently  illximi- 
nating  a  periodic  structure. 

The  Lau  effect  can  form  the  basis  of  an  interferometer  as 

6  7 

first  described  by  Bartelt  and  Jahns  and  later  by  Sudol.  The 
classic  Lau  effect  requires  propagation  from  one  grating  to 
another  to  produce  the  interference  fringes.  The  interferometer 
uses  a  modified  version  where  the  first  grating  is  imaged  onto 
the  second  to  give  an  effective  distance  of  zero.  The  inter¬ 
ference  fringes  are  viewed  in  the  far-field  of  the  second  grating 
by  use  of  a  lens.  The  setup  is  diagrammed  in  Figure  2.  An  inco¬ 
herent  source  (in  our  case  a  tungsten  lamp)  is  Imaged  onto  the 
first  grating  to  provide  an  Incoherent  periodic  source.  is 

imaged  onto  G^  by  the  lenses  L^  and  L^ .  Interference  is  observed 
in  the  focal  plane  of  lens  L^,  which  is  also  the  image  of  plane 
P.  The  interference  phenomenon  can  be  understood  from  a 
geometrical  point  of  view  by  considering  Figure  4.  The  light 
falling  on  G^  will  consist  of  many  plane  waves  at  different 
angles  and  wavelengths.  Any  one  plane  wave  will  be  broken  up 
into  three  major  components  which  we  designate  as  the  0 
(undiffracted),  +1,  and  -1  orders.  These  components  are  Imaged 
onto  G2  where  they  are  further  diffracted.  The  second  diffrac¬ 
tion  causes  certain  orders  to  overlap.  The  +1  order  from  the  -1 
beam,  and  the  -1  order  from  the  +1  beam  both  overlap  with  the 
undiffracted  order.  The  three  orders  come  to  a  common  focus 
after  passing  through  L^ .  It  is  easy  to  see  from  the  diagram 
that  the  system  operates  independently  of  wavelength.  Different 
wavelengths  of  light  will  be  diffracted  different  amounts  by  G^, 
but  the  dispersion  will  be  compensated  at  G^ . 

Light  diffracted  at  G^  comes  to  three  foci  in  the  focal 
plane  of  L^.  If  G^  has  grating  spacing  b,  the  separation  of  the 
foci  is  A  =  XF/b  where  F  is  the  focal  length  of  L^^ .  Because 
these  three  points  originate  from  the  same  beam,  they  are 


coherent  with  respect  to  each  other.  This  is  exactly  the  conclu¬ 
sion  drawn  from  Equation  (7). 

The  Lau  interferometer  is  a  shearing  interferometer;  that 
is,  it  overlaps  adjacent  points  {the  foci  in  Figure  4)  in  the 
object.  As  such  it  displays  the  derivative  of  a  phase  object. 

In  the  absence  of  any  phase  disturbance  between  and  there 
is  no  phase  difference  between  the  adjacent  foci  and  no  inter¬ 
ference  is  observed  in  the  final  plane.  A  quadratic  phase 
disturbance  will  cause  straight-line  fringes  to  appear  in  the 
final  plane.  Straight-line  bias  fringes  can  also  be  formed  by 
displacing  .  If  is  moved  to  some  arbitrary  distance  d  in 
front  of  lens  L^,  the  waves  traveling  between  and  will  no 
longer  be  plane  but  will  have  quadratic  phase  fronts.  The  three 
foci  will  now  have  a  quadratic  phase  relationship  between  them, 
and  bias  fringes  will  appear  in  the  observations  plane.  The 
fringes  will  be  black  and  white  because  the  Interferometer  is 
achromatic. 

Unlike  a  conventional  coherent  Interferometer,  the  fringes 
from  a  Lau  Interferometer  do  not  have  an  Intuitively  obvious 
Interpretation.  The  interpretation  can  be  aided  through  analogy 
with  a  shearing  interferometer  and  by  considering  the  simple  case 
of  bias  fringes  with  no  phase  object  present.  Referring  to 
Figure  2,  grating  G^^  is  moved  a  distance  6  (positive  or  negative) 
away  from  the  front  focal  plane  of  lens  .  No  phase  object  is 
present  at  plane  P.  The  resulting  image  in  the  focal  plane  of 
lens  Lg  is  a  set  of  straight-line  fringes  (see  Figure  6  in 
Section  II. 4).  A  rather  complicated  analysis  of  the  coherence 
function  shows  that  the  intensity  in  the  final  plane  consists  of 
a  shape  function  convolved  with  a  set  of  delta  functions,  which 
is  indeed  a  pattern  of  straight-line  fringes.  For  fringes  to 
occur  the  spacing  b'  of  grating  G^  must  be  related  to  the  spacing 
b  of  grating  G^^  by  b  =  (Fj/F2)b'.  In  other  words,  the  image  of 
grating  must  be  the  same  size  as  grating  G^ . 
condition  is  met  the  fringes  occur  at 


When  this 


where  k  is  an  integer.  The  factor  ^3/^2  magnification 

between  plane  P  and  the  image  plane,  which  we  may  ignore  for  our 
present  purposes. 

We  wish  to  relate  the  spacing  x  in  Equation  8  to  the  phase 
front  in  plane  P.  A  single  point  (which  we  consider  to  be  on- 
axis  for  simplicity)  gives  rise  to  a  wavefront  represented  by 
exp[-i (n/X) (6/F  ^)x  J  in  plane  P.  Because  is  a  grating  there 
is  another  wave  which  originates  at  the  same  point  that  has  a 
linear  phase  shift  proportional  to  the  grating  frequency.  The 
second  wave  has  a  phase  front  at  P  represented  by 
exp[-l(ir/A)  (6/F^)  (x-X(F^/b)  )^J  where  b  is  the  grating  spacing. 
Other  waves  are  also  present  from  the  higher  grating  diffraction 
orders,  but  we  may  ignore  them  here. 

The  grating  overlaps  the  two  waves  so  that  they  add 
coherently  in  the  image  plane.  The  intensity  there  is 
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=  2  +  2cos  [2  X  -  ^  .  (9 

Each  point  on  G^  causes  a  similar  fringe  pattern,  and  all  the 
intensities  are  summed  to  give  the  total  pattern.  The  linear 
phase  terms  can  all  be  ignored.  The  maxima  of  the  pattern  occur 
when 
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where  k  is  an  Integer,  exactly  as  in  Equation  (8). 


We  have  established  that  Lau  interference  patterns  are  ana¬ 
logous  to  shearing  interference  patterns;  therefore,  we  may  use 
the  standard  Interpretation  of  shearing  fringes  for  the  Lau 
interferometer . 

An  interferogram  may  be  written  as  cos  (2?r^(x))  emd  a 
wavefront  as  2iru(x) .  The  shearing  relationship  between  u  and  ^ 


—  S  =  ♦(x) 


where  S  is  the  shear.  For  the  Lau  interferometer 


For  the  present  example  we  have 
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which  is  indeed  the  case  for  a  spherical  wavefront.  This  example 
also  illustrates  why  the  Lau  interferometer  is  achromatic.  A 
phase  front  varies  as  A~^  whereas  the  shear  varies  directly  as  A. 
Therefore,  the  wavelength  cancels  out  of  Equation  (11)  and  <t)(x)  is 
always  independent  of  wavelength. 


4. 


LAU  INTERFEROMETER  EXPERIMENTS 


A  breadboard  Lau  interferometer  was  constructed  to 
demonstrate  its  properties  and  show  its  usefulness  for  wind  tun¬ 
nel  measurements.  The  system  was  identical  to  that  of  Figure  2. 
The  position  of  was  variable,  and  a  phase  object  was  Inserted 
between  lenses  and  L2 ■  All  lenses  had  200  mm  focal  lengths. 
The  gratings  both  had  50  1/in  (2  I /mm)  spacing.  The  light  source 
was  a  tungsten  bulb  in  a  lamp  housing  with  a  variable  aperture 
and  variable  focus  condenser  lens.  The  phase  object  was  a  small 
stream  of  gas  obtained  through  a  hypodermic  needle  attached  to  a 
can  of  freon.  A  very  slight  pressure  on  the  gas  can's  valve  pro¬ 
duced  the  phase  object. 

When  aligned,  and  with  no  phase  disturbance  present,  the 
output  image  has  uniform  intensity.  Close  inspection  shows  the 
image  of  the  needle  to  have  two  ghost  images  associated  with  it. 
These  Images  are  a  result  of  imaging  through  the  second  grating. 
Only  the  +1  orders  from  the  grating  are  evident.  If  the  gratings 
are  moved  slightly  out  of  place,  large  fringes  characteristic  of 
the  shearing  pattern  for  spherical  aberration  can  be  seen.  These 
fringes  have  some  coloration  (red  on  one  edge,  green  on  the 
other).  The  fringes  are  the  result  of  residual  spherical  and 
chromatic  aberrations  in  the  system.  Using  higher-quality  lenses 
would  correct  this  problem. 

Figure  5  shows  an  image  of  the  gas  stream  coming  from  the 
hypodermic  needle.  The  image  is  similar  to  that  derived  from  a 
schlleren  system,  which  is  not  surprising  since  schlleren  is  also 
a  form  of  shearing  interferometry.  The  flow  is  dark  along  the 
bottom  edge  and  bright  along  the  top  edge,  which  is  what  we  would 
expect  for  the  derivative  of  a  flow  with  cylindrical  shape. 

If  Gj  is  moved  away  from  the  front  focal  plane  of  Lj,  bias 
fringes  appear  across  the  image.  Figure  6  shows  how  the  phase 
disturbance  affects  the  bias  fringes.  It  is  easy  to  make  a 


Figure  6. 


Lau  Interfsrogram  of  a  Gas  Stream  When  Bias  Fringes 
are  Added. 


Figure  7. 


Lau  Interferometer  of  Object  Used  in  Figure  6,  but 
Object  Moved  Out  of  the  Region  of  Focus. 


qualitative  evaluation  of  the  phase  flow  just  from  the  pho¬ 
tograph.  It  is  interesting  to  note  that  the  bias  fringes  do  not 
have  uniform  intensity,  but  rather  appear  to  be  darkest  at  the 
edge  with  a  light  stripe  in  the  center.  The  theory  of  the  Lau 
effect  predicts  that  the  interference  fringes  will  have  an  inten¬ 
sity  profile  corresponding  to  the  autocorrelation  of  the  gra¬ 
ting  lines  in  62-  For  our  experiments  G2  was  a  Ronchi  ruling 
with  square  profile  lines,  so  the  autocorrelation  is  a  triangle 
function,  just  as  we  observe.  The  appearance  of  bias  fringes 
from  defocused  gratings  has  never  been  noted  before  in  the 
literature . 

Figure  4  shows  that  for  interference  to  occur  the  gratings 
and  G2  must  be  matched.  More  specifically,  the  periodic 
coherence  formation  at  the  plane  of  the  phase  disturbance  must 
match  the  grating  G2.  We  see  by  Equation  (7)  that  if  Gj  is 
defocused  (afO)  the  period  of  the  coherence  function  changes 
with  the  distance  away  from  the  focal  plane  of  .  If  the  phase 
disturbance  moves  too  far  away  from  the  focal  plane,  the  sampling 
is  no  longer  matched  to  G2  and  interference  no  longer  occurs. 
Therefore,  the  Lau  interferometer  can  be  made  sensitive  to  the 
position  of  the  phase  disturbance.  Figure  7  shows  the  same  phase 
disturbance  as  in  Figure  6  moved  out  of  focus.  It  no  longer  has 
an  effect  on  the  position  of  the  bias  fringes.  Its  only  effect 
is  to  reduce  the  fringe  contrast.  Of  course,  we  might  expect 
some  position-dependent  effect  just  from  the  nature  of  an  inco¬ 
herent  imaging  system.  Indeed,  if  we  increase  the  f /number  of 
the  system  by  stopping  down  lens  Lj ,  the  depth  of  focus  of 
increases  and  the  position-dependent  effect  is  no  longer  so 
apparent.  However,  we  can  definitely  ascribe  some  position- 
dependence  to  the  bias  fringes.  With  no  bias  fringes  across  the 
field  the  phase  object  is  no  longer  visible  after  it  moves  85  mm 
from  the  focal  plane.  With  the  bias  fringes  present  the  object 
is  no  longer  visible  after  a  translation  of  only  35  mm.  So  the 
bias  fringes  have  the  effect  of  narrowing  the  range  over  which  a 
phase  object  is  visible. 


The  defocus  effect  has  an  interesting  implication  for  wind 
tunnel  applications.  Phase  disturbances  will  be  visible  only 
within  a  narrow  region.  Any  other  phase  variations  will  not  be 
visible  as  changes  in  the  bias  fringes.  They  will  only  reduce 
the  overall  fringe  contrast.  However,  a  coherent  Interferometer 
will  Integrate  the  phase  over  the  entire  region  between  Lj  and 
L2 .  Previous  workers  have  only  been  concerned  with  In-f ocus 
images  in  the  Lau  interferometer.  We  are  the  first  to  describe 
the  position-dependent  effects. 

To  formally  establish  the  Interpretation  of  Lau  fringes  we 
used  a  controlled  object  from  which  we  could  get  a  known  phase 
front .  We  wrapped  nlchrome  wire  around  a  glass  tube  to  make  a 
small  heating  filament.  The  heated  wire  is  shown  in  a 
Mach-Zehnder  Interf erogram  in  Figure  8.  The  rod  is  viewed  on  end 
to  show  the  thermal  gradient .  A  Mach-Zehnder  interf erogram  is  a 
direct  measure  of  phase  change,  so  it  is  easy  to  see  how  the 
thermal  gradient  decreases  with  distance  away  from  the  rod. 

The  unheated  rod  is  shown  in  Figure  9  in  the  Lau  inter¬ 
ferometer.  The  bias  fringes  are  aligned  parallel  to  the  rod. 

When  the  rod  is  heated  as  in  Figure  10  the  fringes  shift  to  show 
the  thermal  gradient.  To  emphasize  the  shift  we  have  used  a 
video  processor  to  subtract  the  images  in  the  heated  and  unheated 
states.  The  fringe  shifts  near  the  rod  are  clearly  visible. 

The  wire  Itself  is  also  visible  because  it  expanded  slightly  when 
heated,  but  this  change  does  not  affect  the  fringe  shifts.  The 
fringe  shift  is  large  near  the  rod  and  becomes  smaller  further 
away.  The  same  behavior  of  the  gradient  can  be  deduced  from  the 
Mach-Zehnder  Interf erogram.  An  exact  numerical  comparison  of  the 
interf erograms  would  require  an  expensive  and  time-consuming  ana¬ 
lysis  beyond  the  constraints  of  the  current  contract,  but 
the  qualitative  correspondence  between  the  interf erograms  is 
obvious,  and  we  are  confident  that  further  work  would  establish 
the  numerical  equivalency. 
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Figure  8.  Mach-Zehnder  Interferogram  of  a  Thermal  Gradient 

from  a  Heated  Nlchrome  Wire  Wrapped  Around  a  Glass 
Rod  (~  2  mm  dlam) .  The  Interferogram  was  taken  end 
on.  The  thermal  gradient  Is  clearly  evident. 
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Figure  9. 


Lau  Interferometer  Bias  Fringes  with  Unheated  Wire. 


Figure  10.  Difference  Image  Obtained  by  Electronically 

Subtracting  Lau  Interf erogram  of  Unheated  Object 
from  the  Heated  Object  Interferogram.  The  shift  In 
the  fringe  position  Is  clearly  evident. 


SECTION  III 


COMBINED  LDA/LAU  INTERFEROMETER 
INTERFERANEMOGRAPHY  SYSTEM 

Combining  the  3-D  LDA  system  with  the  Lau  Interferometer 
proceeds  in  a  straightforward  manner.  The  LDA  system  (as  worked 
out  in  Appendix  A)  should  be  formed  from  a  2-D  LDA  system  at 
wavelength  and  a  1-D  LDA  system  at  X^  combined  at  +45o  and 
-45<>  to  the  normal  to  the  principal  flowfield  velocity  direction. 
The  Lau  interferometer  can  be  set  up  in  line  along  the  normal 
(Figure  11)  giving  each  optical  system  clear  access. 

To  provide  real-time  data  acquisition  the  phase 
shift  in  the  Lau  interferometer  observation  plane  S  (Figure  11) 
will  need  to  be  measured  electronically.  The  detector  must  be 
positioned  in  the  plane  S  at  the  location  where  the  LDA  volume  is 
focused.  To  interpret  the  signal  unambiguously,  it  is  necessary 
that  the  actual  signal  intensity  be  measured,  i.e.,  DC  coupled. 
This  can  lead  to  problems  resulting  from  a  variation  in  the  light 
source  intensity  and/or  changes  in  optical  transmission. 
Therefore,  a  reference  detector  should  be  used  to  monitor  all 
background  fluctuations,  and  the  resulting  signals  should  be 
ratloed . 

Another  problem  is  the  27t  phase  ambiguity.  This  can  be 
overcome  by  dithering  the  position  of  along  the  plane  of  the 
grating.  The  fringes  will  vibrate  and  the  relative  phase  of  the 
vibrator  and  the  AC  component  of  the  detector  signal  will  deter¬ 
mine  the  direction  of  the  detected  phase  change  and  whether  or 
not  a  27r  phase  transition  occurred. 

The  integrated  LDA/Lau  inter feranemography  system  will  be 
able  to  provide  real-time  data  acquisition.  Subsequent  computer 
interpretation  will  be  required  for  both  signals.  As  an  added 
benefit,  at  the  expense  of  more  data  channels,  additional  detec¬ 
tors  can  be  placed  in  the  Lau  detection  plane,  and  density  data 
of  a  more  global  nature  will  be  available  to  assist  interpreta¬ 
tion  of  the  Interference  pattern. 


SECTION  IV 

SUMMARY  AND  CONCLUSIONS 

The  Lau  interferometer  Is  a  true  coherent  interference 
device  (as  opposed  to  moire)  which  needs  only  incoherent  illumi¬ 
nation.  It  has  an  in-line  configuration  which  uses  no  beam 
splitters  or  spatial  filters.  It  is  also  extremely  stable. 
(Striking  the  table  with  a  rubber  mallet  did  not  produce  any 
fluctuation  in  the  fringes.) 

The  Lau  is  a  shearing  interferometer,  so  it  presents  the 
derivative  of  the  phase  variation.  It  can  be  used  in  either  a 
schlleren  mode  or  with  bias  fringes  to  aid  in  visualizing  the 
phase  disturbance.  With  a  low  f /number  system  using  bias 
fringes,  the  region  of  phase  visualization  is  limited  to  the 
focal  region  of  the  first  lens.  The  Lau  interferometer  is 
capable  of  visualizing  one  plane  in  a  thick  phase  disturbance. 

We  must  emphasize  that  our  understanding  of  the  Lau  inter¬ 
ferometer  to  this  point  is  largely  qualitative.  To  fully 
understand  and  apply  it  an  in-depth  theoretical  and  experimental 
study  will  be  required.  Several  questions  must  be  answered  by 
this  study.  The  primary  objective  is  interpreting  the  fringe 
patterns.  They  contain  derivative  information,  but  the  infor¬ 
mation  is  not  in  a  simple  form.  In  the  configuration  used  in 
this  study  the  output  consists  of  three  overlapped  images,  rather 
than  two  as  in  other  shearing  interferometers.  It  is  also 
necessary  to  determine  under  what  conditions  a  simple  fringe 
interpretation  is  valid.  We  anticipate  that  a  weak  phase  distur¬ 
bance  is  required,  but  that  the  conditions  are  influenced  by  the 
f /number  of  the  system  and  the  frequency  of  .  Likewise,  these 
variables  affect  the  depth  of  focus  of  the  phase  imaging. 

Along  with  the  theoretical  aspects  of  the  interference  pat¬ 
tern  the  engineering  aspects  of  the  entire  system  must  also  be 
examined.  The  choice  of  gratings  is  critical  to  the  system. 


What  happens  if  radial  gratings  are  used,  rather  than  linear  gra- 

O 

tings?  Previous  work  has  shown  that  in  classic  moire  inter¬ 
ferometry  a  phase  grating  may  be  used  as  the  second  grating.  If 
a  thick  phase  grating  is  used  for  G^,  it  might  eliminate  one  of 
the  overlapping  images,  thus  simplifying  the  overall  image.  Also 
a  large-aperture,  low  f /number  system  should  be  set  up,  perhaps 
using  reflective  optics,  to  provide  a  system  better  suited  for 
wind  tunnel  work.  We  feel  that  a  research  and  development 
program  to  address  these  issues  will  result  in  a  valuable  new 
tool  for  wind-tunnel  diagnostics. 

The  LDA  system  is  well  in  hand.  The  system  is  thoroughly 
understood  and  the  technology  required  to  use  LDA  is  commercially 
available.  The  shortcoming  with  some  commercial  applications  of 
3-D  systems  is  in  the  on-axls  velocity  resolution.  By  using  the 
wider  angle  system,  discussed  in  Appendix  A,  the  system  geometry 
can  be  tailored  to  provide  the  desired  performance.  The  short¬ 
comings,  where  there  are  any,  tend  to  be  in  the  data 
acqulsition/slgrnal  processing  systems  and  in  optical  access. 
Currently,  only  customized  computer  systems  have  the  capability 
to  provide  the  high-speed  data-gathering  power  required  for  scale 
length  and  correlation  time  measurements. 
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SECTION  A-1 
INTRODUCTION 


As  part  of  Contract  No.  F33615-83-C-3020  the  University  of 
Dayton  has  conducted  a  study  of  technologies  involved  and 
required  for  an  inter feranemography  system.  The  current  tech¬ 
nologies  of  interferometry  and  laser  Doppler  anemometry  were 
evaluated  for  their  applicability  to  interf eranemographic  tech¬ 
niques  as  well  as  other  techniques  such  as  Rayleigh  scattering 
and  holographic  mapping  of  flow-field  particles. 

The  specified  requirements  for  an  interferanemographlc 
system  are  to  provide  real-tlme/on-line  simultaneous  measurements 
of  density  and  velocity  in  complex  flow  fields. 

A  generally  accepted  method  for  measurement  of  air  density 
is  interferometry.  Within  the  past  decade  holographic  inter¬ 
ferometric  techniques  have  been  used  successfully  to  obtain  den¬ 
sity  data  in  a  variety  of  wind  tunnel  studies.  However,  the 
general  application  and  wide  use  of  interferometry  is  hindered  in 
wind  tunnel  studies  by  the  following:  (a)  single  llne-of-sight 
views  of  the  flow  field  limit  the  application  to  2-D  or  3-D 
axisymmetric  fields;  (b)  insufficient  measurements  for  optically 
weak  fields  (l.e.,  hypersonic  flows  or  low  density  flows);  (c) 
manual  data  reduction  methods  now  widely  used. 

Laser  Doppler  techniques  have  been  used  to  measure  veloci¬ 
ties  in  subsonic  and  supersonic  flow  fields.  Laser  Doppler 
Anemometry  (LDA)  is  a  process  whereby  changes  in  the  frequency 
(Doppler  shifting)  of  a  laser  radiation  scattered  from  particles 
moving  through  an  optical  control  volume  are  monitored  to  compute 
the  velocity  of  the  moving  particles.  The  particles  are  kept 
small  (0.2  to  1.0  micron  in  diameter),  typically  so  that  their 
velocities  are  the  same  as  the  local  velocity  of  the  fluid 
stream. 


Section  A-2  of  this  report  presents  a  development  of  the 
requirements  and  techniques  for  interferometric  measurements  of  a 
general  three-dimensional  flow  field.  As  a  part  of  the  study 
using  interferometry,  it  was  necessary  to  characterize  the  inter¬ 
ferometer  sensitivities  for  the  different  wind  tunnels  at 
Wright-Patterson  Air  Force  Base,  and  the  results  of  this  com¬ 
parison  are  presented  in  Section  A-2. 

Section  A-3  presents  a  development  of  the  Laser  Doppler 
Anemometry  requirements  and  techniques  for  characterization  of 
local  velocities  in  a  complex  flow  field.  Alternate  techniques 
for  measurement  of  the  local  velocities  are  compared  to  LOA  to 
see  if  they  offer  any  advantage  over  LDA  in  work  by  the  Flight 
Dynamics  Laboratory. 

Section  A-4  summarizes  the  potential  of  existing  technology 
to  obtain  real-tlme/on-llne  interferanemographic  measurements  in 
complex  flow  fields. 


SECTION  A-2 

INTERFEROMETRY  MEASUREMENT  OF  FLOW-FIELD  DENSITIES 

Aerodynamic  measurements  of  density  using  Interferometry 
require  the  use  of  optical  equipment  amd  processes  which  are 
capable  of  recording  and  Interfering  light  waves  of  special 
character.  Several  processes  Involving  holography  have  been 
developed  and  used  successfully  In  wind  tunnel  applications,  and 
these  results  are  documented  In  the  literature.  In  the  present 
study,  the  prospects  for  obtaining  real-time  Interferometric  data 
simultaneously  with  laser  veloclmetrlc  measurements  are  sought. 
Motivation  for  the  development  of  this  capability  comes  from  atn 
aerodynamic  need  to  be  able  to  measure  the  perturbation  factors 
comprising  the  Reynolds  shear  stresses  In  air  flows.  These  fac¬ 
tors  are  the  perturbations  In  both  the  density  and  velocity. 

• 

An  Important  physical  condition  to  understand  In  the  use  of 
Interferometry  Is  the  fact  that  the  measurable  fringe  shifts  of 
an  Inter ferogram  result  from  a  summation  of  all  changes  In  den¬ 
sity  along  the  path  traversed  by  the  light  waves.  This  means  the 
values  of  density  calculated  from  Interferometric  measurements 
are  path-dependent  quantities.  Furthermore,  unless  the  flow 
field  Is  two-dimensional  and  the  optical  cucls  of  the  Inter¬ 
ferometer  Is  aligned  with  the  axis  of  the  flow  along  which  the 
fluid  properties  are  Invariant,  more  than  one  view  of  the  flow 
field  Is  required  to  obtain  density  data.  It  Is  regrettable  that 
this  limitation  places  definite  restrictions  on  the  general  use¬ 
fulness  of  Interferometry  in  wind  tunnel  applications,  because 
many  facilities  in  use  today  simply  prohibit  multiple  viewing  of 
the  f low . 

In  response  to  the  present  study,  a  general  description  of 
the  interferometer  process  applicable  to  wind  tunnel  studies  is 
presented  along  with  a  survey  on  Interferometric  sensitivities 
for  the  specific  wind  tunnel  facilities  of  the  Air  Force  Wright 
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Aeronautical  Laboratories,  Flight  Dynamics  Laboratory  (AFWAL/FI). 
Also  given  Is  a  description  of  a  net*  Interferometric  technique. 
Point  Diffraction  Interferometry  (PDI),  which  appears  to  be  a 
promising  approach  to  obtaining  "real  time”  Interferometric  data. 


A-2.1  BASIC  INTERFEROMETRY  FOR  FLOW  FIELDS 

The  Interference  fringes  obtained  with  Interferometry  are 
related  to  the  refractive  Index  of  the  gases  In  the  test  section 
and  the  reference  section  by 

S(x,y,z)  =  ^  /  [n<x,y,z)-n^^j]dl  (A-1) 

where  S(x,y,z)  Is  the  measurable  fringe  shift,  X  Is  the  wave¬ 
length  of  the  light  and  n(x,y ,z)-nj.ef  is  the  relative  change  In 
refractive  Index. 

For  gases,  the  Index  of  fraction  Is  related  to  density 
through  the  Dale  Gladstone  constant  Kqq  as, 

^  n(x,y,z)-l 
DG  p(x,y,z) 

Therefore,  Equation  (A-1)  becomes 

S(x,y,z)  =  G  /  f(x,y,z)  dl 

1 

where 


Is  a  constant ,  and 

f(x,y,z)  =  _  1  (A-5) 

^ref 

Equation  (A-3)  Is  difficult  to  resolve  because  the  Integrand, 
f(x,y,z),  depends  on  three  spatial  coordinates  (see  Figure  A-2.1) 
When  applied  to  a  plane,  x  -  x^  where  Xq  denotes  a  constant,  the 
dependence  Is  reduced  to  two  variables  and  Equation  (A-3)  becomes 


(A-2) 

(A-3) 


Figure  A-2.1  Illustration  of  Interferometric  view  of  a  general 
3-D  field. 


Figure  A-2.2  Ray  orientation  for  2-D  axlsymmetrlc  field. 


S(y,z)  =  G  J  f  (y,z)  dl 

I 


(A-6) 


No  further  general  simplification  Is  possible  except  for  cases 
where  the  flow  field  Is  either  2-D  or  3-D  axlsymmetrlc .  For 
these  cases 


S(y)  =  G  [f(y)]L 


(2-D) 


(A-7) 


S(y)  =  G  / 

y 


'  2  2 
r  -y 


(3-D  axlsymmetrlc)  (A-8) 
(See  Figure  A-2.2) 


The  density  Is  obtained  by  solving  these  relationships  for  the 
density  function.  This  results  In 


-  GL^  *  "ref 


(2-D) 


(A-9) 


4  ^ 


dy  +  p. 


r  [y2_^2 


(3-D  axlsymmetrlc).  (A-10) 


For  asymmetric  fields,  a  general  solution  to  the  basic 
equation  does  not  exist.  Instead,  multiple  views  of  the  field 
are  required  In  conjunction  with  the  use  of  special  transfor¬ 
mations.  For  steady  flows,  rotation  of  the  model  can  be  used  to 
obtain  multiple  views,  thereby  permitting  computation  of  the 
density  from  the  optical  data.  However,  for  unsteady  flows,  the 
model  rotation  technique  is  not  truly  accurate,  and  results  In 
considerable  uncertainty  In  the  data.  Exact  solutions  are 
possible  when  the  field  Is  viewed  over  180  degrees  In  increments 
of  approximately  5  to  8  degrees.  Smaller  viewing  ranges  can  be 
used  (generally  not  less  than  90  degrees) ,  but  error  Is  Inherent 
In  these  results. 


A-2.2  INTERFEROMETRIC  SENSITIVITIES  FOR  AFWAL/FI  HIND  TUNNELS 


As  a  means  of  evaluating  Interferometric  measuring  capabil¬ 
ities  in  the  AFWAL/FI  wind  tunnels,  the  Interferometric  sen¬ 
sitivity  is  used  to  estimate  representative  fringe  shift 
magnitudes .  These  results  are  presented  in  Table  A-2 . 1 .  Since  the 
length  factor  in  the  interferometric  integral  equation  varies 
from  one  application  to  another,  an  interferometric  sensitivity 
per  unit  length  is, 
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f(y)L  X 

For  the  present  evaluation,  the  reference  density  is 

P  P  — ^ 

■  -sS-  N1  +  M  ^ 
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assuming  Isentroplc  flow  from  the  wind  tunnel  nozzle  to  the  test 
rhombus.  The  zero  subscripts  denote  the  total  state,  and  R,  y 
and  Mqo  are,  respectively,  287  J/Kg-*K,  the  ratio  of  specific 
heats  «  1.4,  the  free  stream  mach  number;  Kdq  And  X  are  the  same 
values  as  previously  reported.  For  wind  tunnel  facilities  of 
reasonably  constant  Tq  and  fixed  Mqo,  variation  of  the  total 
pressure  is  seen  to  directly  affect  the  interferometric  sen¬ 
sitivity;  the  order  of  magnitude  spread  for  the  Mach  3  facility 
illustrates  this  effect. 


To  better  compare  the  interferometric  measuring  capabilities 
of  these  facilities,  a  15-cm  length  is  used  to  obtain  the  actual 
interferometric  sensitivities.  These  results  occur,  for  example, 
if  a  15-ca  wide  flat  plate  model  were  tested  in  each  of  the 
facilities.  Finally,  the  most  Important  result,  the  actual 
fringe  shifts,  are  shown  in  the  far-right  column  of  Table  2.1. 
Note  that  the  results  for  the  4- in.  HWT  are  parenthesized, 
because  a  IS-cm  model  can  not  be  tested  in  this  facility. 


The  fringe  shifts  presented  in  Table  A-2.1  are  based  upon  the 
magnitude  of  the  density  change  across  a  two-dimensional  tur- 
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bulent  boundary  layer  on  a  flat  plate.  The  span  of  the  plate  Is 
considered  to  be  15  cm,  and  the  plate  is  assumed  to  be  at  zero 
angle  of  attack  in  the  uniform  steady  floM  air  streeuB  of  the  test 
rhombus.  For  this  situation,  the  density  function  is 


f(y) 


1  - 
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1  - 


oo 
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since  the  static  pressure  variation  through  the  boundary  layer  is 
constant.  The  static  wall  temperature,  varies  with  time  as 

the  model  remains  in  the  flow.  Depending  on  the  heat  transfer 
rates,  the  wall  condition  could  reach  an  adiabatic  wall  state,  at 
which  time 
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where  r^.  is  the  recovery  factor. 
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for  turbulent  air  flow.  Therefore,  assuming  the  adiabatic  wall 
condition  exists,  the  density  function  becomes 


f(y)  -  1  - 
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The  usefulness  of  these  results  is  apparent  in  Table  A-2 . 1 .  Both 
high  Reynolds  number  facilities  (Mach  3  and  6  wind  tunnels) 
appear  good  for  Interferometric  studies,  the  TGF  is  marginally 
good  for  transonic  and  supersonic  studies,  and  the  hypersonic 
facilities  are  essentially  no  good  for  Interferometric  measure¬ 
ments.  The  fringe  shifts  (49.73  and  33.76)  of  the  20-ln.  HHT 
and  the  9- in.  SAW  are  too  high  and  would  probably  result  in  a 
general  lack  of  fringe  distinction. 


As  previously  mentioned,  the  field  width,  L,  is  an  important 
factor  in  computing  interferometric  sensitivity  and  the 


corresponding  fringe  shifts,  and  the  15-cm  length  of  these  results 
serves  only  as  a  basis  for  comparison.  One  notices  that  as  L 
decreases,  the  fringe  shift  diminishes,  and  In  fact  it  becomes 
small  even  for  the  high  density  flows  of  the  Mach  3  and  6 
facilities  when  L  is  on  the  order  of  a  few  millimeters.  This 
might  be  true  for  boundary  layer  studies  on  cone  or  cylinder 
models. 


A-2.3  INTERFEROMETRY  FOR  REAL-TIME  MEASUREMENT  OF  DENSITIES 

Present  work  In  interferometry  at  the  Flight  Dynamics 
Laboratory  is  not  capable  of  providing  real-time  density  measure¬ 
ments  for  flow  analysis  in  the  study  of  dynamic  events.  At  this 
time  the  interference  patterns  are  recorded  during  a  test  and 
then  evaluated  in  a  non-real-time  mode  using  mamual  data  reduc¬ 
tion  techniques. 

As  shown  in  Section  A-2.1,  if  interferometry  is  used  to 
evaluate  density  profiles  in  complex  asymmetric  flows,  multiple 
views  of  the  field  must  be  recorded.  In  order  to  record  and 
evaluate  in  real-time  density  profiles  using  interferometry,  an 
instrument  is  required  that  includes  the  following  components: 

1.  Interferometer  Illumination  Source, 

2.  Interferometric  Device, 

3.  Interference  Pattern  Recording  System, 

4.  Interference  Pattern  Evaluation  System, 

5.  Algorithms  for  Determination  of  Densities, 

It  was  originally  proposed  that  holographic  interferometric 
techniques  would  be  evaluated  for  its  feasibility  to  be 
integrated  into  an  inter feranemographlc  system. 

A-2.3. 1  Holographic  Interferometry 

A  hologram  is  a  photographic  recording  of  a  coherent  object 
wavefront  and  a  reference  wavefront.  When  a  processed  hologram 


Is  relllumlnated  by  the  reference  wavefront,  the  hologram  will 
diffract  the  light  In  the  reference  beam  and  reproduce  an  exact 
replica  of  the  object  wavefront.  After  reconstruction,  the 
reconstructed  wavefront  of  the  flow  field  can  be  used  for  study 
of  the  flow  field.  By  using  a  pulsed  laser  source  to  record  the 
hologram,  the  record  freezes  In  time  a  picture  of  the  flow  field 
as  It  existed  during  the  laser  exposure  time.  With  laser  pulses 
of  50  ns  or  less  the  temporal  variations  of  the  flow  field  can  be 
observed.  By  comparing  a  time  history  of  the  Interference 
between  the  object  wavefront  and  a  reconstruction  of  the  object 
wavefront  as  It  existed  In  some  reference  state,  quantitative 
Information  can  be  obtained  of  the  temporal  variations  In  the  gas 
densities  of  a  flow  field.  Figure  A-2.3  shows  a  schematic  drawing 
of  a  candidate  holographic  Interferometer  system  for  use  In 
studying  flow  fields.  Optical  windows  are  used  to  provide  optical 
access  to  the  test  region.  For  real-time  evaluation  of  the 
Interference  pattern  between  a  reconstructed  reference  object 
wavefront  and  the  object  wavefront  during  a  test,  a  television 
camera  would  be  used  to  view  the  real-time  Interference  pattern. 
The  laser  source  would  be  pulsed  In  sequence  with  the  30-Hz  tele¬ 
vision  frame  rate,  and  the  digitized  Interference  pattern  from 
the  television  camera  would  be  evaluated  to  obtain  the  gas  den¬ 
sities  as  explained  In  Section  A-2.1. 

A-2.3. 1.1  Sources 

Although  a  pulsed  laser  would  not  be  required  to  record  the 
reference  object  wavefront,  It  will  be  necessary  for  real-time 
evaluation  during  a  test  run.  The  pulse  length  of  the  laser  must 
be  short  enough  to  freeze  the  temporal  effects  and  have  good 
enough  coherence  and  beam  quality  to  maintain  a  well-behaved 
coherent  wavefront  through  the  flow  field  volume  relative  to  the 
reference  beam.  To  obtain  good  historical  Information  of  a  flow 
field  In  a  dynamic  flow  requires  a  laser  source  with  the  short 
pulses  of  a  Q-swltched  ruby  laser  and  the  versatility  and  con¬ 
tinuous,  adjustable  repetition  rate  capabilities  of  a  modulated 


argon  laser.  The  two  most  likely  sources  available  are  an 
exclmer  laser  pumped  dye  and  a  metal  vapor  laser.  A  standalone 
exclmer  laser  Is  not  a  usable  source  because  of  the  very  high 
transverse  mode  character  of  the  beam.  Multiple  transverse  modes 
will  tend  to  scramble  the  coherence  of  the  beam  across  the  field, 
making  holography  difficult  If  not  Impossible. 

Exclmer  lasers  emit  only  In  the  ultraviolet,  have  no  spa¬ 
tial  coherence,  and  are  limited  to  500  Hz,  although  the  pulse 
energies  are  high  (100  to  200  mJ)  and  the  pulse  length  Is  short 
(20  nsec).  The  exclmer  laser  could  be  converted  to  a  usable 
light  source  by  using  It  to  pump  a  dye  laser.  By  pumping  a  dye 
laser  with  an  exclmer  laser,  the  short,  hlgh-repetltlon  pulse  of 
light  of  the  exclmer  can  be  maintained  while  potentially 
Increasing  the  coherence  of  the  beam  through  the  properties  of 
the  dye.  The  dye  laser  will  emit  at  many  wavelengths  through  the 
visible  spectrum  and  can  be  restricted  to  give  a  good  coherence 
length.  The  pulse  energy  with  an  exclmer  pumped  dye  Is  an  order 
of  magnitude  less  than  the  exclmer  alone  (10  to  30  mJ)  and  the 
pulse  length  Is  still  very  short. 

Because  It  is  a  new  laser  system,  there  Is  not  very  much 
experimental  data  on  the  beam  quality  and  coherence  produced  by 
combining  the  unacceptable  exclmer  beam  with  the  well-behaved  dye 
laser  medium.  Pulse  repetition  rates  are  also  somewhat  limited 
at  this  time  (about  500  Hz)  although  higher  rates  with  sufficient 
energy  will  probably  be  available  In  the  future.  The  cost  of 
these  systems  Is  currently  high  ($60  K  to  $60  K) ,  because  of  the 
two  lasers  required.  The  physical  size  of  these  systems  Is  also 
quite  large,  requiring  an  area  about  4  m  by  1  m  by  1  m.  If  lower 
repetition  rates  are  acceptable,  Molectron  Corporation  produces  a 
lower  cost,  compact  design  system. 

There  are  a  number  of  other  possible  lower  cost  lasers  that 
provide  lower  repetition  rates.  A  dye  laser  alone  has  a  much 
longer  pulse  length  than  an  exclmer  (up  to  500  nsec)  but  Is 


limited  in  repetition  rates  to  under  50  Hz.  If  high  energies 
(over  100  mJ)  are  required,  dye  lasers  are  limited  to  under  10 
Hz.  A  dye  laser  is  adequate  when  viewing  with  a  video  system 
(which  limits  the  repetition  rate  to  30  Hz  because  of  the  video 
frame  rate),  and  they  are  low  in  cost  (typically  under  $20  K  from 
Candella  Corp.).  Frequency-doubled,  electro-optically  Q-swltched 
Nd:YAG  will  put  out  the  same  energy  as  the  dye  laser  in  a  much 
shorter  pulse  (15  nsec)  but  it  is  currently  limited  to  30  Hz 
repetition  rates  and  it  costs  twice  as  much  as  a  dye  laser.  A 
program  is  underway  with  NASA/Lewis  attempting  to  push  this  repe¬ 
tition  rate  up  for  the  Nd:YAG  lasers  so  that  it  is  suitable  for 
combustion  diagnostic  applications. 

The  Alexandrite  laser  can  produce  energies  similar  to  those 
of  the  dye  laser  (100  mJ)  with  a  shorter  pulse  length  but  at  only 
slightly  higher  repetition  rates  (40  Hz).  Because  of  the  higher 
cost  and  limited  coherence  properties  of  the  Alexandrite  laser  it 
is  not  as  attractive  an  option  as  the  dye  laser. 

Frequency-doubled,  acousto-optically  Q-swltched  NdrYAG 
will  operate  at  repetition  rates  up  to  50  kHz  with  a  200  nsec 
pulse  length,  but  at  pulse  energies  of  less  than  5  mJ.  There  has 
been  very  little  work  done  with  frequency-doubled,  Nd:YAG 
acousto-optically  Q-swltched,  so  the  stability  and  coherence  pro¬ 
perties  are  not  well  characterized.  It  is  very  likely  that  with 
development  work  a  higher  energy,  better  controlled  frequency- 
doubled,  acousto-optically  Q-switched  Nd:YAG  could  be  built.  The 
laser  would  need  to  use  amplifier  systems  and  feedback  circuitry 
which  would  make  it  fairly  expensive.  NdzYAG  technology  is  well 
established  though,  and  it  is  a  reasonably  reliable  system.  The 
primary  limitation  of  the  acousto-optically  Q-switched  Nd:YAG  would 
be  its  long  pulse  length  which  would  limit  its  application  in 
holographic  recording. 

The  repetition  rate  capability  of  the  laser  system  will 
directly  affect  the  measurement  and  mapping  accuracy  of  a  par¬ 
ticle  system. 


A  metal  vapor  laser  system  Is  reliable  and  self-contained, 
and  It  has  greater  versatility  than  the  exclmer  systems.  In 
addition  to  Its  having  a  good  working  coherence  length,  by  using 
an  externally  heated  system  the  University  of  Dayton  has  shown^ 
that  the  metal  vapor  laser  can  be  prompted  to  pulse  at  low  rates 
of  a  few  Hz  to  around  5  kHz  with  the  potential  of  operating  up  to 
50  kHz  In  short  bursts.  The  synchronizing  electronics  for  such  a 
system  Is  similar  to  that  used  for  the  strobed  argon  laser  devel¬ 
oped  by  the  University  of  Dayton  Research  Institute  (UDRI).^ 

The  metal  vapor  laser  Is  a  system  with  properties  similar 
to  those  of  the  acousto-optlcally  Q-swltched,  frequency-doubled 
Nd:YAG  In  terms  of  repetition  rate  and  scalability  of  pulse 
energy  but  with  a  much  shorter  (25  nsec)  pulse  length  and  better 
Inherent  coherence  and  stability  properties.  We  evaluated  the 
manganese  metal  vapor  laser  with  a  UDRI-bullt  breadboard  laser 
shown  In  Figure  A-2.4.  Initial  studies  with  this  laser  Indicate 
that  It  can  be  built  to  operate  at  1  Hz  up  to  a  few  hxindred  Hz 
continuous  and  operated  In  long  bursts  of  perhaps  100  pulses  up 
to  10  kHz  and  potentially  up  to  50  kHz.  The  pulse  length  Is 
still  very  short  and  the  coherence  length  Is  greater  than  5  cm. 
The  breadboard  laser  put  out  over  100  |iJ  per  pulse.  A  larger 
copper  laser  can  output  up  to  12  mJ  with  current  designs  (larger 
systems  are  also  possible).  For  a  self-heated  metal  vapor  laser, 
the  metal  temperature  Is  not  dependent  on  a  high  repetition  rate 
to  sustain  operation.  Therefore,  the  laser  can  be  cued  to 
operate  from  single  pulse  to  a  few  hundred  or  a  thousand  Hz. 

This  Is  currently  the  only  laser  of  this  type  available  with  a 
pulse  length  on  the  order  of  20  to  40  nsec.  The  pulse  energy, 
coherence,  and  stability  exhibited  by  this  system  are  superior  to 
those  of  any  comparable  laser  system  at  this  time. 

As  with  many  newer  types  of  lasers  (dyes,  exclmers,  etc.), 
the  metal  vapor  laser  Is  a  hlgh-malntenance  Item;  the  laser  must 
be  replaced  or  recycled  regularly.  The  warmup  time  required  has 
also  been  viewed  as  a  drawback.  Originally,  a  long  warmup  time 


was  necessary  to  reduce  the  chance  of  tube  breakage.  A  new  tube 
design  may  eliminate  this  requirement  (which  is  currently 
required  by  most  commercial  metal  vapor  lasers) .  The  warmup  time 
of  the  breadboard  laser  was  found  to  be  1-1/2  to  2  hours. 

Due  to  the  low  average  power  of  the  breadboard  system 
(about  10  mw,  primarily  unpolarlzed)  this  particular  system  may 
not  be  practical  as  a  source  for  making  holograms  of  large 
volumes.  However,  even  this  breadboard  laser  produces  more 
energy  per  pulse  than  an  argon  laser  at  the  pulse  length  which 
seemed  desirable  to  freeze  dynamic  vibrations  of  a  few  hundred 
Hz  (typically  100  to  200  )sec  with  the  argon  laser).  It  will 
be  possible  to  make  a  hologram  of  a  field  using  a  metal  vapor 
laser,  then  view  the  reconstruction  with  a  small  argon  laser 
(which  is  very  long-lived  and  has  low  maintenance) . 

Table  A-2.2  reviews  the  potential  lasers  for  holographic  flow 
studies  based  upon  discussions  with  manufacturers  and  University 
of  Dayton  experience.  The  pulse  energy  required  will  depend  on 
the  type  of  film  used  to  make  the  hologram.  For  typical  expo¬ 
sures  this  is  likely  to  be  less  than  1  mJ.  To  record  the 
hologram  it  is  also  necessary  to  have  a  coherence  length  suf¬ 
ficient  to  cover  all  path  length  variations  such  as  the  marginal 
rays  going  through  the  edge  of  the  collection  lens  or  due  to  the 
tilt  of  the  reference  beam.  For  practical  application,  that 
would  imply  a  coherence  length  of  5  cm  or  greater  (a  linewidth  of 
0.0035  nm  or  less).  From  Table  A-2.2  it  is  seen  that  several 
existing  lasers  are  suitable,  and  the  actual  laser  selected  would 
be  determined  by  the  budget  and  the  measurements  to  be  made. 

A-2.3.1.2  Recording  Material 

An  important  factor  which  will  affect  the  quality  of  a 
hologram  and  dictate  the  laser  pulse  energy  required  is  the  film 
used  to  record  the  hologram.  There  are  a  number  of  different 
emulsions  available  on  film  and  plates,  notably  from  Kodak  and 


viewing 


Agf a-Gavaert ,  which  are  capable  of  recording  a  hologram.  These 
films  vary  in  sensitivity  and,  inversely,  in  resolution.  These 
films  are  summarized  in  Table  A-2.3.  The  speed  of  any  of  these 
films  can  be  increased  by  a  factor  of  3  to  10  by  preprocessing 
the  film,  either  by  pre-exposlng  it  with  a  clean  uniform  beam  or 
by  predeveloping  (except  for  the  Laser  Technology  film  which  is 
already  preprocessed) .  This  preprocessing  typically  will  not 
affect  the  quality  of  the  hologram  (since  the  entire  emulsion  is 
affected  tiniformly)  but  may  reduce  the  efficiency  of  the  hologram 
by  reducing  the  contrast  of  the  recording.  A  lower-contrast 
hologram  will  be  less  noisy  than  a  high-contrast  hologram  since  a 
high-contrast  hologram  will  accentuate  every  light  variation, 
like  turning  up  the  contrast  on  a  video  system. 

The  Laser  Technology  film  is  effectively  already  prepro¬ 
cessed  since  it  is  exposed  wet.  When  we  spoke  with  Laser 
Technology  regarding  this  application,  they  mentioned  that  their 
holograms  are  typically  fairly  noisy  due  to  the  very  active 
developer  they  use.  Their  process  is  oriented  toward  fast  deve¬ 
loping  of  holograms  primarily  for  field  applications. 

A-2.3. 1.3  Interference  Pattern  Recording  and  Computer  Evaluation 

To  obtain  real-time  measurements  of  the  flow-field  density 
distribution  a  television  camera  would  observe  the  real-time 
interference  pattern,  and  the  digitized  video  image  would  be  com¬ 
puter  analyzed  to  determine  the  fringe  pattern.  The  interference 
fringe  pattern  Information  would  then  be  used  to  calculate  the 
refractive  index  and  density  variations  through  the  flow  field. 

Software  requirements  for  processing  the  image  of  a  fringe 
pattern  start  with  the  utility  routines  needed  to  control  the 
image  digitizer  and  then  transfer  the  digitized  image  into  memory 
and  finally,  if  desired,  to  the  display  screen.  Other  utility 
routines  are  needed  to  enhance  the  contrast,  zoom  in  on  a 
selected  section  of  the  image,  mask  out  defects,  and  perform 
other  elementary  image-processing  functions. 


TABLE  A-2.3 


HOLOGRAPHIC  RECORDING  MATERIALS 


Sensitivity 

in  Green/cm^  (Gamma) 


Material 

(Density  =  1) 

Resolution 

Contrast 

Kodak 

649.  F 

plate 

800 

ergs 

2000 

+ 

5 

SO-343 

film 

1000 

ergs 

2000 

+ 

5 

120-02 

plate 

500 

ergs 

2000 

+ 

5 

SO-173 

film 

500 

ergs 

2000 

+ 

4 

Mlnlcard  II 

SO-424 

film 

50 

ergs 

1250 

4 

131-101 

plate 

30 

ergs 

1250 

7 

SO-253 

film 

30 

ergs 

1250 

7 

Aafa-Gavaert 

8E56HD-AHI 

plate/film 

350 

ergs 

2500 

+ 

4 

8E75HD 

plate/film 

800 

ergs 

2500 

+ 

3 

10E56 

plate/film 

20 

ergs 

1500 

7 

10E75 

plate/film 

120 

ergs 

1500 

4 

Laser  Technology 


Wet  processing  film 
Laser  test  film 


8  ergs 
(red  tested) 


2000  + 


lo  in  lo  '4' 


The  next  task  is  to  match  the  image  to  the  important 
features  of  the  flow  field  area  and  test  part.  The  coordinate 
system  of  the  digitized  image  must  be  made  to  correspond  to  the 
actual  object  space.  This  is  done  by  manually  Identifying  key 
features,  such  as  corners  of  objects  or  external  fiducial  marks, 
whose  location  in  object  space  is  knotm.  Software  will  be  pro¬ 
vided  to  outline  the  shape  of  any  test  fixtures  in  the  image. 

This  gives  the  operator  positive  assurance  of  the  accuracy  of  the 
coordinate  transformation  so  that  he  may  halt  the  analysis  pro¬ 
cess  at  that  point  if  the  display  does  not  look  right. 

An  interactive  setup  dialogue  is  the  best  approach  for  the 
learning  phase  of  the  system  to  define  the  shapes  of  test  fix¬ 
tures  and  models.  The  result  is  a  data  file  that  can  be  used  to 
identify  and  measure  flow  fields  around  the  model  in  a  test.  An 
operator  may  be  required  to  locate  these  features  on  each  image. 
Simple  geometric  models  can  be  defined,  so  that  this  iden¬ 
tification  procedure  will  require  only  a  few  seconds  of  work 
with  a  cursor. 

Locating  fringe  centers  Is  difficult,  €uid  no  one  algorithm 
is  suited  for  all  types  of  fringe  patterns.  Over  uniform  or 
background  regions  of  the  image,  the  fringes  may  be  straight  and 
evenly  spaced.  A  coarse  scan  across  the  image,  with  large  steps 
between  scans,  may  provide  enough  information  to  describe  the 
model  geometry  in  these  regions  by  a  low-order  polynomial 
approximation.  These  regions  may  be  digitized  rapidly.  In  some 
regions  of  the  image  the  fringe  pattern  may  be  complex  and  low  in 
contrast.  Such  regions  require  high-resolution  but  relatively 
slow  fringe-following  techniques.  In  highly  complex  areas  of  the 
image,  the  software  can  enlarge  the  region  of  interest  and  label 
locations  of  known  fringes.  For  high-resolution  measurement 
applications,  an  operator  may  be  required  to  untangle  the  pattern 
manually  with  the  use  of  a  joystick  or  digitizing  tablet.  The 
computer  can  assist  by  modifying  the  display  in  response  to 
operator  commands. 


Once  fringe  centers  have  been  located,  a  final  data  file  Is 
prepared  and  saved  for  analysis.  Absolute  model  measurement 
requires  that  full  complexed  fringe  data  be  obtained  and  analyzed 
for  all  areas  of  Interest.  If  a  difference  fringe  pattern  Is 
created  to  show  only  deviations  from  the  master  part,  this  auialy- 
sls  can  be  made  much  simpler.  It  Is  thought  that  the  time 
required  for  evaluation  of  a  single  Interference  pattern  would  be 
about  2  minutes.  The  cost  and  details  of  this  type  of  system 
will  be  presented  to  the  Air  Force  when  desired.  A  conceptual 
diagrauH  of  a  developmental  fringe  analysis  system  is  shown  in 
Figure  A-2 . 5 . 

A-2.3.2  Point  Diffraction  Interferometry 

A  relatively  new  method  for  obtaining  interferometric  data 
of  flow  fields  involves  the  use  of  Point  Diffraction  Interferometry 
(PDI).  First  reported  by  Smartt  in  1975, ^  PDI  offers  the  use  of 
a  single-pass,  relatively  vibration** free  interferometer  that 
could  be  used  for  continuous  surveillance  of  the  changes  in  the 
optical  path  lengths  of  light  waves  traversing  the  flow.  Unlike 
interferometers  of  the  usual  type,  PDI  is  easier  to  use,  because 
only  one  light  beam  is  required  to  make  the  measurements,  and 
while  not  yet  done,  application  to  wind-tunnel  testing  using  con¬ 
ventional  schlieren  optical  systems  seems  appropriate. 

In  the  conceptual  use  of  a  wind  tunnel  PDI ,  a  beam  of 
coherent  light  passes  through  the  test  section  of  the  tunnel 
using  an  optical  scheme  comparable  to  a  Toepler  schlieren  system. 
Once  the  beam  is  captured  by  the  second  schlieren  mirror  (or 
lens),  the  light  waves  are  directed  into  the  PDI.  As  illustrated 
in  Figure  A-2. 6,  the  optical  components  of  the  PDI  resemble  those 
of  a  conventional  Mach  Zehnder  interferometer.  One  set  of  light 
waves  passes  through  one  path  of  the  interferometer  unaltered 
while  the  other  set  is  focused  and  filtered  similarly  to  the  way 
coherent  radiation  can  be  spatially  filtered  so  as  to  transmit 
selected  bemds  of  diffracted  radiation.  Here,  the  filtering  is 
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done  so  as  to  pass  only  the  rays  which  are  not  refracted  by  the 
density  field  in  the  wind  tunnel.  As  in  the  conventional  Mach 
Zehnder  Interferometer,  the  beams  are  reunited  by  the  second 
splitter  plate  which  is  also  used  to  regulate  the  structure  of 
the  interference.  From  here,  the  usual  fringe  reading  and  appli¬ 
cable  data  reduction  processes  can  be  used  to  compute  flow-field 
densities  from  the  optical  data. 

The  prospects  of  PDI  are  excellent  because  of  its  vibration 
insensitivity.  As  with  holographic  interferometry  a  pulsed 
source  and  television  camera  could  be  used  for  recording.  Thus 
far,  however,  PDI  has  not  been  applied  to  wind  tunnel  tests. 

A-2.3.3  Heterodyne  Interferometry 

Both  the  holographic  interferometer  and  the  point  diffrac¬ 
tion  interferometer  can  be  operated  as  heterodyne  Interferometers 
to  provide  real-time  Interf erometry . ^  In  the  heterodyne  mode  a 
frequency  shifter  (l.e.,  moving  mirror  in  the  reference  beam)  is 
used  to  modulate  the  phase  in  the  reference  beam  in  a  periodic 
fashion  to  provide  real-time  Interferometry.  Since  heterodyne 
interferometry  provides  a  direct  measurement  of  the  object  wave 
phase  differences,  it  eliminates  the  need  to  find  and  evaluate 
fringe  pattern  positions  as  described  in  the  interference  pattern 
evaluation  system.  Section  A-2.3.1.3.  This  is  especially  impor¬ 
tant  in  the  analysis  of  flow  fields  where  the  fringe  pattern  can 
be  badly  distorted  near  the  edges  of  test  objects  in  the  flow 
fields.  Likewise  since  it  is  a  direct  measurement  technique  the 
amount  of  computer  processing  is  much  less,  since  the  only  require¬ 
ment  is  to  compute  the  densities  and  refractive  indices  from  the 
directly  measured  phase  variations  in  the  object  wavefront.  With 
this  type  of  instrument  the  phase  measurements  are  instantly 
obtained  instead  of  the  2  minutes  in  the  system  described  in 
Section  A-2.3.1.3  where  the  fringe  patterns  are  digitized.  A 
detector  array  with  the  required  spatial  resolution  is  used  to 
directly  measure  the  real-time  phase  differences.  It  is  expected 


that  It  would  cost  at  least  as  much  as  the  systems  of  Section 
A-2.3.1.3  because  no  work  has  been  done  with  this  technique  for 
flow-field  evaluation. 


SECTION  A-3 

FLOW  FIELD  VELOCITY  MEASUREMENTS 


Optical  techniques  have  been  extensively  utilized  for 
nonintrusive  probing  of  fluid  velocities.  A  variety  of  different 
approaches  have  been  used  with  the  major  efforts  concentrated  in 
either  time-of-f light  systems  or  in  measurement  of  the  Doppler 
shift  of  the  frequency  of  scattered  light.  The  time-of-f light 
techniques  can  measure  the  time  to  cross  a  single  light  beam 
( one-spot ),S  whose  cross-section  is  known,  or  they  can  measure 
the  time  to  cross  between  two  light  beams  ( two-spot ),6  whose 
spacing  is  known.  Tlme-of-f light  techniques  are  characterized  by 
strong  signals,  the  light  beams  are  tightly  focused,  and  the 
techniques  can  make  measurements  where  signal  interferences  such 
as  background  luminosity  or  surface  scattering  would  stop  other 
techniques.  Their  drawback  is  that  they  can  only  measure  one 
part  of  the  velocity,  either  a  single  component  or  the  magnitude 
of  the  velocity  in  a  plane.  Extension  of  this  method  to  multi- 
component  measurement  is  not  practical. 

Doppler  shift  techniques  for  measuring  velocity  have  a  long 
history,  especially  in  the  field  of  astronomy  where  the  Doppler 
shift  is  used  to  measure  the  recession  speed  of  stellar  bodies. 
These  speeds  are  very  high  and  the  spectral  shifts,  due  to  the 
Doppler  effect,  are  large  and  easily  discernible,  even  with 
coarse  instrumentation.  As  the  speed  decreases,  the  Doppler 
shift  decreases  proportionally  and  much  more  sophisticated 
instrumentation  is  required  to  accurately  measure  the  change  in 
frequency.  Commonly  Fabry-Perot  Interferometers  (FPI)  are  used  for 
high-resolution  measurements.  Laser-based  Doppler  measurement 
instruments  have  been  constructed,  using  FPI ' s  either  at  the  source 
or  at  the  detection  end. 

The  use  of  a  laser  as  a  source  of  radiation  to  interrogate 
the  moving  body  has  eliminated  the  need  to  have  a  body  that 
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emitted  radiation,  as  in  the  stellar  observations.  The  laser's 
bandwidth  can  be  much  narrower  than  the  half widths  associated 
with  spectral  transitions,  and  this  feature  has  been  used  in 
several  techniques  to  provide  measurements  of  quite  small  veloci¬ 
ties,  useful  for  laboratory  aerodynamic  studies.  One  approach 
has  been  to  scatter  the  laser  radiation  from  a  solid  body  moving 
with  the  fluid.  The  scattering  is  elastic  and  the  scattered 
radiation  will  be  Doppler  shifted  in  frequency.  If  the  body  is 
moving  fast  enough,  so  that  the  shift  is  greater  than  the  laser's 
linewidth,  the  FPI  can  be  used  to  measure  the  velocity  field. ? 
Another  technique,  used  in  wind  tiuinels,  has  been  to  seed  the 
flow  with  Na  atoms.*  The  Na  can  absorb  and  re-emlt  the  laser 
radiation,  the  "scattering"  strength  proportional  to  the  Na 
absorption  coefficient  at  the  laser  wavelength.  Irradiating  the 
flow  at  an  angle  to  the  main  flow  direction  will  result  in  a 
shift  of  the  absorption  band,  according  to  the  Doppler  shift, 
from  where  it  would  be  in  a  stationary  sample.  By  sweeping  the 
frequency  of  a  laser's  narrow  linewidth,  and  Illuminating  both  a 
stationary  sample  and  the  moving  fluid,  the  Doppler  shift  can  be 
determined  directly  from  the  intensities  scattered  from  the  two 
samples. 8  Zimmerman  and  Mlles^  have  developed  a  clever  modifica¬ 
tion  that  uses  a  retroref lector  and  eliminates  the  need  for  the 
stationary  sample. 

*  This  is  a  difficult  technical  feat  and  special  care  must  be 
taken  to  prevent  the  formation  of  dimers.  Also,  this  technique 
is  only  well  suited  to  helium  flowflelds  since  He  does  not  quench 
excited  Na. 
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Another  approach  is  to  seed  the  flowfield  with  iodine.  The 
I2  spectrum  is  very  rich  in  the  visible  region  and  several  lines 
overlap  Argon  laser-lines.  By  using  an  etalon  (effectively  an 
FPI  in  the  source  instead  of  at  the  detector),  the  laser  can  be 
made  to  oscillate  on  a  single  longitudinal  mode  and  this  mode  can 
be  swept  over  the  laser  bandwidth.  As  a  result,  the  I2  absorp¬ 
tion,  and  subsequent  fluorescence  intensity,  will  depend 
upon  the  location  of  the  laser  line  within  the  Doppler  shifted 
linewidth  of  the  I 2  absorption  line,  and  this  I 2  line  can  be 
mapped  by  shifting  the  laserlO"^^  (the  frequency  chemges  in 
Integral  amounts  of  a  fixed  quantity  which  can  be  accurately 
determined  independently) .  An  alternative  approach^^^  has  been  to 
observe  the  scattering  from  iodine  using  four  beams,  the  forward 
and  retro-ref lected  beams  in  two  directions  to  the  flow.  The 
authors  attempted  to  determine  two  components  of  the  flow  velo¬ 
city  from  the  measurements.  In  any  technique  requiring  multiple 
observations, » 13®  the  observations  must  all  occur  within  one 
correlation  time  of  the  flow  in  order  that  the  observations  can 
be  interpreted  as  a  flow  velocity  measurement.  Unfortunately  the 
time  restriction  in  most  aerodynamic  flows  of  interest  is  too 
short  for  the  two-dimensional  mapping  probes  to  handle  and,  as  a 
result,  with  current  technology  these  probes  cannot  provide 
meaningful  measurements  of  velocity  fluctuation. 

The  last  Doppler  technique  that  will  be  introduced  is  the 
measurement  of  the  Doppler  shift  of  elastically  scattered 
radiation  using  either  heterodyne  or  homodyne  mixing. This  is 
the  most  commonly  employed  nonlntrusive  velocity  measurement 
technique,  and  instrumentation  is  available  from  several  commer¬ 
cial  sources.  The  principle  of  the  technique  makes  use  of  the 
temporal  coherence  (monochromaticity)  of  the  laser  source.  The 
Doppler-shifted  scattered-radiation  frequency  is  down-shifted  to 
where  it  can  be  measured  electronically.  The  down-shift  fre¬ 
quency  is  usually  derived  from  the  same  laser  source  and  the 
Doppler  shift  can  be  measured  precisely,  enabling  an  accurate 
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determination  of  the  velocity  in  the  flowfield.  The  major  limi¬ 
tation  is  the  frequency  range  that  can  be  measured  so  that  the 
method  becomes  less  useful  at  higher  Mach  numbers.  The  technique 
requires  elastic  scattering,  usually  from  solid  bodies  in  the 
flow,  and  is  known  as  either  laser  Doppler  anemometry  (LDA)  or 
laser  Doppler  velocimetry  (LDV)  depending  on  how  accurately  the 
solids  are  following  the  fluid  flow.  LDA  techniques  can  be 
extended  to  multidimensional  measurements  in  a  straightforward 
manner  and  therefore  it  is  the  favored  technique  in  aerodynamic 
velocity  measurements.  Consequently,  in  the  following  sections 
the  LDA  technique  and  its  application  to  boundary  layer  measure¬ 
ment  with  a  supersonic  free  stream  flow  will  be  studied  in 
detail.  This  is  the  most  trying  condition  for  boundary  layer 
measurements  and  the  system  designed  for  this  flow  will  also  work 
well  in  lower-speed  tunnels. 

In  the  sections  that  follow  a  detailed  design  and  error 
analysis  of  a  3-D  LDA  will  be  worked  out  based  upon  the  dual-beam 
LDA.  This  type  of  system  represents  the  current  state-of-the-art 
multi-component  velocity  measurement  instrument  used  in  wind 
tunnels.  The  criteria  for  selecting  the  configuration  will  be 
presented  in  a  general  format.  The  detailed  examples  will  be 
relevant  to  the  Flight  Dynamics  Laboratory's  LDA  system  in  opera¬ 
tion  in  the  transonic  wind  tunnel. 

A- 3.1  LASER  DOPPLER  ANEMOMETRY 

In  the  far  field,  where  the  scattered  light  wave  is  a 
spherical  wave,  the  scattered  spherical  field  can  be  described  in 
terms  of  its  frequency 

=  1^0  i  ^  (A-12) 

where  Pg  is  the  scattered  field's  frequency,  Uq  'the  Incident 
frequency,  X  is  the  incident  wavelength  in  the  medium,  v  is  the 


vector  velocity  of  the  scatterer,  and  f  and  §  are  unit  vectors  in 
the  direction  of  observation  and  the  incident  laser  direction 
respectively  (Figure  A-3.1).  The  Doppler  shift,  (Vs-Vq),  depends 
linearly  on  the  component  of  v  in  the  direction  (f-S) .  It  is 
Important  to  be  aware  of  the  relative  sizes  of  the  parameters. 
Typically  Vg  ~  ^  10^^  Hz  when  A  =  514.5  nm  (Argon  green 

line),  so  that  the  maximum  Doppler  shift  would  be  -o  3.9 
MHz/(m/s).  Even  at  high  velocities,  the  Doppler  shift  is  a 
small  fraction  of  the  Incident/scattered  frequency  making  direct 
observation  difficult  (ex.  3.9  GHz  9  1000  m/s  represents  a  shift 
of  6.7  parts  per  million  of  the  original  frequency). 

A-3.1.1  Downbeat inq 

It  is  possible  to  beat  the  scattered  radiation  with  another 
light  beam  derived  from  the  same  laser  source.  The  mixing  takes 
place  in  the  photodetector,  a  device  that  responds  to  intensity, 
not  to  amplitude.  In  a  simplified  example  consider  two  sinusoidal 
fields  sin  U]^t  and  sin  U2t  Incident  on  the  detector.  The  detec¬ 
tor  would  respond  to  [sin  w^t  +  sin  (i)2t)2.  By  use  of  trigono¬ 
metric  identities  this  becomes  (sin  «it  +  sin  «2t]2  =  i  + 

(cos  2wit  +  cos  2u2t)  +  cos(«i-«2)'*^  “  cos(wi  +  W2)t.  The  fre¬ 
quencies  2u>i ,  2<m)2,  and  ui  +  «i>2  large  that  the  detector 

cannot  respond  and  it  averages  the  result  to  zero.  If  (i>i-(i)2  is 
sufficiently  small  the  detector  can  follow  this  component's 
variation.  The  photodetector  downmixes  the  signal  to  the  fre¬ 
quency  {ui-u)2)  • 


There  are  three  optical  configurations  used  in  LDA  to  pro¬ 
vide  a  reference  frequency. These  are  known  as  the  reference 
beam  mode,  the  dual  scatter  mode,  and  the  dual-beam  mode,  all 
shown  in  Figure  A-3.2.  In  the  reference  beam  mode  a  weak  reference 
beam  is  split  from  the  main  beam  and  mixed  with  the  scattered 
radiation  at  the  detector.  The  disadvantage  to  this  mode 
(Figure  A-3.2. a)  is  that  the  wider  the  collection  aperture  (to 
Increase  signal),  the  greater  the  uncertainty  in  r  and  there- 


(a)  Reference  beam  mode. 


(b)  Dual  scatter  node. 
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(c)  Dual  beam  mode. 


Figure  A-3.2  LDA  optical  configurations. 


fore  the  greater  the  luicertalnty  in  the  velocity  measurement. 

Also,  the  alignment  requirements  for  heterodyning  are  severe  and 
difficult  to  maintain. The  advantage  is  that  the  reference 
beam  can  be  brought  from  any  angle  so  that  any  direction  of  velo¬ 
city  measurement  can  be  selected. 

The  dual-scatter  mode  mixes  together  the  light  scattered  in 
two  different  directions  (Figure  A-3.2.b).  The  resulting  Doppler 
signal  yj)  =  Vg!  ~  •^s2  v«(ri-r2)  is  Independent  of  the  direction 

of  the  incident  beam.  The  disadvantage  of  this  mode  is,  as  with 
the  reference  beam  mode,  that  the  aperture  must  be  restricted  in 

order  to  define  rj  and  r2  with  sufficient  accuracy  that  the  com¬ 
ponent  of  9  in  the  direction  (^i-^2)  can  be  computed  to  the 
desired  accuracy.  This  is  a  severe  restriction  on  the  collection 
aperture. 

The  dual-beam  mode  mixes  the  scattered  radiation  from  two 
input  beams,  whose  unit  propagation  vectors  s^  and  ^2'  that 
cross  at  an  angle  0.  The  resulting  scattered  frequencies  from 
each  beam,  in  a  particular  direction  r,  are 
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and 
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(A-13b) 


As  a  result  of  the  mixing,  the  Doppler  signal  generated  by  the 
detector  Is 

(A-14) 
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which  is  independent  of  r.  This  is  the  great  advantage  of  the 
dual  beam  method;  the  Doppler  frequency  resulting  from  mixing  the 
scattered  radiations  is  the  same  everywhere  in  space.  Con¬ 
sequently  the  scattered  radiation  can  be  collected  over  a  large 
aperture,  and  a  substantially  stronger  signal  will  result  than  for 
either  of  the  other  LDA  methods.  The  Doppler  frequency  is 
linearly  proportional  to  the  component  of  velocity  in  the  direc¬ 
tion  S-  -  S,  with  the  relation 


where 


(A-16) 


K  =  ^  sin  e/2 

A-3.1.2  "Fringe"  Model 

When  two  plane  parallel  laser  beams  with  the  same  polariza¬ 
tion  Intersect  on  a  screen,  a  set  of  Interference  fringes  occurs 
with  the  spacing 
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Consequently  researchers  proposed  the  "fringe"  model  to  explain 
the  LOA  signal.  The  Doppler  frequency  results  from  the  scat¬ 
tering  of  a  particle  as  it  crosses  a  series  of  interference 
fringes.  This  model  provides  a  simple  conceptual  picture  to 
explain  the  LDA  process  which  is  correct  only  in  limited 
situations  and  should  not  be  used  in  a  general  sense. But  so 
long  as  the  scatterer  size  is  substantially  less  than  the 
"fringe"  spacing  the  results  will  be  accurate. 


A  difficulty  arises  with  the  use  of  Equation  (15)  for  genera¬ 
lized  turbulent  flows.  The  Doppler  signal  is  not  direction  sen¬ 
sitive,  i.e.,  the  detector  interprets  and  -t-vi)  as  the  same 
frequency.  This  is  overcome  by  a  single-sideband  modification  in 
which  the  frequency  of  one  of  the  incident  beams  is  shifted  by  a 
fixed  amount,  n3c,  usually  using  an  acousto-optic  Bragg  cell. 

The  Doppler  signal  now  becomes 

^D  ®  *^BC  ^  Kv  (A-18) 

and  so  long  as  the  negative  value  of  Kv  never  exceeds  n^c  the 
Doppler  signal  can  be  interpreted  unambiguously.  Difficulties 
arise  in  selecting  an  optimal  value  of  ngc  ^or  measurements 
through  flows  that  have  large  mean  velocity  changes,  such  as 
botindary  layers  in  supersonic  wind  tunnels.  Often  a  compromise 
value  must  be  selected  to  provide  adequate  LDA  performance  in  all 


regions.  An  example  illustrating  a  selection  of  V3Q  is  worked 
through  later. 


Using  the  concept  of  the  "fringe”  model,  several  operational 
requirements  for  the  optical  system  become  evident.  The  Doppler 
signal  will  be  generated  by  scattering  from  the  region  in  space 
where  the  two  incident  beams  cross.  For  the  "fringes"  to  have 
the  same  spacing  throughout  this  region  it  is  necessary  that  the 
beams  be  plane  waves.  Also,  for  ease  of  Interpreting  the  flow 
Information,  the  region  of  space  in  which  the  beams  cross  should 
be  very  small.  Fortunately  these  two  requirements  can  be 
fulfilled  simultaneously  by  focussing  the  two  input  beams  cuid 
crossing  them  at  their  waists.  The  resulting  beam  parameters 
(e~2  intensity  profiles,  wavefront  radius  as  a  function  of  posi¬ 
tion,  etc.)  have  been  presented  in  the  literature^^ > and  will 
not  be  developed  here.  Only  the  results  will  be  used. 

A-3.1.3  Dual  Beam  LDA 

A  simple  configuration,  which  will  provide  the  required 
optical  arrangement,  is  shown  in  Figure  A-3.3.  The  key  element  is 
the  lens  which  both  crosses  the  beams  and  focusses  them  at  the 
crossing  point.  In  actual  practice  a  telescope  should  be  used  to 
ensure  that  the  beams  focus  at  the  crossing  location  and  the 
interference  pattern,  as  seen  on  a  screen,  should  be  examined 
throughout  the  crossing  region  (the  LDA  spot).  The  beams  will 
focus  to  a  diameter 


d  _2  =  (4/7T)  A  (f/D  )  (A-19) 
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where  f  is  the  focal  length  of  the  lens,  D  _o  is  the  diameter  of 

e 

the  beams  at  the  lens,  and  f/D  o  is  the  f-number  of  the  focused 
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beam.  The  LDA  spot  is  an  ellipsoid  with  dimensions 
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and  the  number  of  fringes  within  the  spot  is 
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It  should  be  noted  that  when  a  Bragg  cell  is  used  the  concept  of 
the  number  of  fringes  has  little  utility.  Instead  the  spot 
crossing  time  and  the  Doppler  frequency  must  be  computed  and 
their  product  will  give  the  maximum  number  of  Doppler  cycles  in 
the  signal . 


A-3.1.4  Multicomponent  Velocity  Systems 

There  are  three  common  methods  for  extending  the  LDA  system 
to  two-dimensional  measurements  and  there  are  several 
configurations  of  each.  The  further  extension  to  three- 
dimensional  measurements  is  more  restrictive.  The  three  two- 
dimensional  dual-beam  LDA  systems  use  either  color,  polarization, 
or  frequency  to  differentiate  the  components.  The  systems  can  be 
either  three-beam  or  four-beam  systems  and  are  shown  as  a  matrix 
in  Figure  A-3.4.  The  optical  system  can  be  rotated  about  the 
centerline  to  provide  a  measurement  coordinate  orientation  better 
suited  to  the  flow  system  being  studied,  or  to  provide  better 
access  around  models.  If  the  beam  pattern  is  translated  with 
respect  to  the  centerline  or  if  the  beams  are  tilted  relative  to 
the  third  axis.  Information  on  the  velocity  component  in  this 
direction  may  be  mixed  into  the  signals  (see  Appendix) .  In  a  two 
dimensional  system  care  must  be  exercised  to  ensure  that  the 
maximum  contribution  of  this  component  to  the  signals  is  less 
than  the  uncertainty  in  the  measurements  themselves. 


The  color  separation  and  polarization  separation  schemes 
have  been  used  extensively.  The  frequency  separation  scheme 
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suffers  from  the  fact  that  the  frequency  range  for  each  component 
is  severely  limited  since  the  two  components  must  be  separated  by 
frequency  filters.  The  extension  from  a  two-dimensional  system 
to  a  three-dimensional  system  usually  entails  the  addition  of  a 
separate  one-dimensional  system  to  a  two-dimensional  system. 

When  appropriate  optical  access  is  available  the  third  component 
can  be  aligned  orthogonally  to  the  other  two  ouid  the  measurement 
is  straightforward.  Typical  systems  use  Argon  ion  lasers  and 
will  combine  polarization  and  color  separation  to  get  all  three 
components  using  the  two  strongest  Argon  lines  (514.5  nm  and 
488.0  nm)  or  they  will  use  three  different  Argon  lines  (514.5  nm, 
488.0  nm  and  either  496.5  nm  or  476.5  nm) .  The  two  strong  lines 
have  comparable  powers  but  the  next  strongest  lines  are  reduced 
to  f  15%  of  the  stronger  lines'  power.  As  a  result  measuring  the 
third  component  in  a  three-color  system  requires  added  care. 

More  often  than  not  the  flow  system  being  measured  does  not 
provide  optical  access  permitting  orthogonal  measurements.  There 
are  two  generic  schemes  in  use  based  upon  measurements  in  a 
nonorthogonal  coordinate  frame.  The  commercial  realizations  of 
these  systems  are  shown  in  Figure  A-3.5.20,21  jn  addition  there  are 
several  individual  systems,  similar  in  appearance  and  operation 
but  using  different  separation  schemes . 22 , 23  one  feature  of  the 
system  shown  in  Figure  A-3.5(a)  is  that  the  scattered  signals  are 
collected  on  the  opposite  arms  from  the  incident  beams.  This 
ensures  that  the  two  arms  are  aligned  to  the  same  spot  and  that 
the  measured  signal  on  each  channel,  at  one  instant  in  time, 
comes  from  the  same  scatterer.  This  is  of  particular  importance 
when  measurements  are  taken  in  a  nonorthogonal  coordinate  system 
or  when  cross-product  data  (e.g.,  Reynolds  stress)  are  desired. 


A-3.1.5  Considerations  in  Practical  Application 


There  are  several  aspects  of  the  acquisition  and  interpreta¬ 
tion  of  LDA  data  that  need  to  be  taken  into  account  in  the  design 
and  operation  of  the  LDA.  First,  there  are  the  uncertainties 


Figure  A-3.5a  Schematic  of  DISA  3-D  LDA  system. 


Figure  A-3.5b  Schematic  of  TSI  3-D  LDA  system  composed  of  a 

3-beam  2-D  frequency  separation  system,  measuring 
two  non-orthogonal  velocity  components,  and  a 
2-beam  1-D  LDA  measuring  the  velocity  component 
orthogonal  to  the  plane  of  the  other  two  com¬ 
ponents  . 
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that  result  when  the  data  are  measured  In  one  coordinate  frame 
and  then  transformed  to  a  different  frame.  Then  there  is  the 
question  of  how  many  measurements  are  required  to  characterize 
the  distribution  to  any  specified  moment.  Both  of  these 
questions  are  relevant  to  aerodynamic  measurements  in  general  auid 
to  the  AFWAL/FI  Transonic  Wind  T\mnel  studies  in  particular. 

A-3. 1.5.1  Coordinate  Transformation 


Consider  a  flowfield  in  which  the  coordinate  frame  of 
Interest  is  x,y,z,  an  orthogonal  set.  Suppose  an  LDA  system  is 
used  that  makes  measurements  along  the  coordinates  x',y',z',  as 
shown  in  Figure  A-3. 6.  The  transformation  from  the  measurement 
frame  to  the  reference  frame  is  given  by: 
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and  is  {-l)^'*'^  multiplied  by  the  determinant  of  the  submatrix 

of  D  formed  by  striking  out  the  row  and  column  containing  the 
(i,j)  element.  For  example: 
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Figure 


Orthogoridl  reference  coordinate  frame  (X,Y,Z)  and 
LDA  coordinate  frame  {X',Y',Z'). 


Z' 


Figure  A-3.7  Relationship  of  measurment  frame  and  reference 
frame . 
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Consider  the  geometry  shown  In  Figure  A-3.5(a).  The  coordinate  frame 
is  such  that  z  and  z'  coincide  and  x,y,x' ,y'  are  as  shown  in 
Figure  A-3.7.  The  transformation  becomes 
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The  uncertainties  associated  with  each  velocity  calculated  (one 
Standard  deviation)  can  be  expressed  as 


«V..  = 
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In  most  systems  of  practical  interest  the  precision  in  deter¬ 
mining  <f>  is  sufficient  that  first  term  [of  order  (fi*)')^]  is 
negligible  no  matter  how  small  the  turbulence  intensity.  If  the 
flowfield  were  to  move  the  beams  about  so  that  this  motion  became 
the  determining  uncertainty  in  0,  then  the  relative  importance  of 
the  terms  might  change.  The  uncertainties  67^  and  6Vy  will  be 


considered  as  resulting  entirely  from  £Vx<  and  6Vyi .  For  the 
geometry  shown  in  Figure  A-3.7, 


6Vx  =1.05  RMS 

(6V) 

(A-30a) 

«Vy  =3.52  RMS 

(6V) 

(A-30b) 

where 


RMS  (6V)  =  -§  /  .  (A-31) 

Therefore,  the  uncertainty  in  the  determination  of  Vy  is  3.5 
times  larger  than  the  RMS  uncertainty.  For  the  geometry  shown  in 
Figure  A-3.5(b)  the  factor  increases  to  6.4x.  It  is  evident  from 
this  discussion  that  the  uncertainty  in  the  velocity  deter¬ 
mination  must  be  reduced  to  as  little  as  practical  and,  also,  in 
a  nonorthogonal  reference  frame,  the  angles  between  the  axes 
should  be  as  close  to  right  angles  as  is  practical.  In  some 
experimental  configurations  the  beams  are  canted  so  that  z  and  z' 
do  not  align  either.  This  further  complicates  the  issue  as  the 
V2 >  component  is  mixed  into  and  Vy.  Analysis  of  this 
situation  is  a  direct  extension  of  the  material  presented.  It 
should  be  emphasized  that  the  uncertainty  in  each  component  is  a 
function  of  the  magnitude  of  the  RMS  uncertainty.  Consequently, 
even  if  the  flow  is  principally  in  one  direction,  the  uncertainty 
in  the  velocity  measurement  of  the  spanwise  component  may  be 
greatly  Influenced  by  the  uncertainty  in  the  streamwise  component. 

A-3 .1.5.2  Counter-Processor  Errors 

The  uncertainties  expressed  in  Equations  (A-29a-b)  include 
only  the  random  uncertainties  resulting  from  measurement  limita¬ 
tions.  Bias  errors  are  not  included.  It  is  presumed  that  careful 
experimentalists  will  ensure  that  the  contribution  resulting 
from  potential  bias  error  sources  is  smaller  than  the  random 
fluctuations.  A  source  of  the  random  fluctuations  will  now  be 
considered  and  the  magnitude  will  be  computed  for  a  "typical" 


LDA  measurement  system.  The  parameters  used  will  be  for  applica¬ 
tion  to  measurements  through  the  boundary  layer  and  into  the  free 
stream  in  a  supersonic  wind  tunnel.  Assumed  is  a  system  with  a 
counter-processor  (C-P) ,  1  ns  resolution,  determining  the  Doppler 
frequency.  The  laser  source  mode  spacing  is  taken  to  be  140  MHz 
so  that  the  maximum  Doppler  frequency  will  be  70  MHz.  (An  etalon 
can  be  used  to  eliminate  the  other  longitudinal  modes  but  70  MHz 
will  remain  a  practical  limit  due  to  the  response  of  the  typical 
photodetector .  For  high-power  Argon  lasers  the  longitudinal  mode 
spacing  is  much  less.)  The  size  of  the  LDA  spot  diameter  will  be 
taken  to  be  250  )m  (d,g)  and  the  free  stream  velocity  will  be  1000 
m/s.  Because  the  crossing  time 

tc  =  dm/v  =  250  ns  (A-32) 

is  minimal  at  the  highest  velocity,  the  Doppler  frequency  must  be 
the  highest  possible  for  the  system  at  this  velocity,  70  MHz. 

The  maximum  number  of  Doppler  cycles  scattered  will  be 


Nc  =  yp  ^c  *  1®  (A-33) 

necessitating  the  C-P  analysis  circuitry  to  be  set  for  2®  cycles 
if  it  is  desired  to  make  efficient  use  of  the  scattering  events, 
most  of  which  will  scatter  less  than  cycles.  If  the  uncer¬ 
tainty  in  the  period  determination  is  +1  C-P  clock  cycle,  then 
5^0  is  0.875S<(  and 
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Will  depend  on  the  frequency  of  the  Bragg  cell.  A  typical  Bragg 
cell  frequency  is  40  MHz.  If  this  were  used  then  for  Uj)  to  be  70 
MHz  at  1000  m/s  would  require  K  =  30  kHz/(m/s),  equivalent  to  a 
"fringe"  spacing  of  33.3  4m.  Therefore 
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at  X  =  514.5  nm,  and  £V  %  2%.  The  spot  length,  Ig),  is  too  long 
to  define  the  scattering  event  location  precisely.  Also,  the  velo 
city  uncertainty  deteriorates  rapidly  as  the  velocity  decreases 
(6V  »  7.7*  when  V  =  100  m/s),  and  the  uncertainty  in  the  orthogo¬ 
nal  component  is  even  larger.  The  situation  can  be  improved  by 
reducing  ngQ  to  as  small  a  value  as  practical.  For  exaunple.  If  a 
dual  Bragg  cell  arrangement  is  employed  to  effecv  1*30  =  5  MHz* 
then  df  becomes  15.4  (im  and  l|g  becomes  3.8  mm.  Again  6);3  = 

0.875*  but  now  6V  =  0.94*.  When  the  velocity  drops  to  100  m/s, 
tj;  =  2.5  |is,  Uj)  =  11.5  MHz,  and  Nc  =  29.  Leaving  the 

C-P  set  for  2^  cycles  results  in  a  frequency  uncertainty  6V3  of 
0.14*  and  6V  =  0.25*.  Even  at  v  =  10  m/s,  61*3  =  0.07*  and 
5V  =  0.6*.  All  the  errors  computed  should  be  acceptable  in  most 
research  situations.  The  system  deteriorates  rapidly  when  two- 
dimensional  measurements  are  required  and  one  measurement  vector 
is  orthogonal  to  the  streamwise  component.  In  this  case  the 
spanwise  velocity  is  zero,  the  Doppler  frequency  is  1^30  but  the 
crossing  time  is  determined  by  the  streamwise  velocity,  yielding 


=  1  at  V  =  1000  m/s,  a  totally  unacceptable  result. 


Either  v. 


can  be  raised  back  to  40  MHz  with  the  accompanying  increase  in 
the  uncertainties,  or  the  optical  system  can  be  rotated  to  +45® 
about  the  streamwise  component.  In  the  latter  case,  the  Doppler 
frequency  will  decrease  if  the  "fringe"  spacing  is  maintained  or 
the  "fringe"  spacing  can  be  reduced  to  maintain  the  Doppler  fre- 


*In  boundary  layer  flow  measurements  V30  can  be  made  small  since 
the  Bragg  cell  offset  frequency  is  only  used  to  ensure  sufficient 
Doppler  cycles  for  analysis  during  the  passage  of  the  scatterer. 
In  recirculating  flows,  especially  where  the  reverse  velocity  is 
comparable  to  the  free  stream  velocity,  V30  must  be  maintained 
near  the  middle  of  the  available  Doppler  bandwidth. 


quency  at  70  MHz  (preferable) .  When  the  system  is  further 
expanded  to  three-dimensional  measurements  using  either  of  the 
schemes  shown  in  Figures  A-3.5(a)  or  A-3.5(b),  the  1-D  system's  optimum 
alignment  will  be  as  parallel  to  the  free  stream  as  possible. 

The  uncertainties  at  lower  velocities  can  be  reduced  by 
increasing  the  number  of  C-P  cycles  used  in  the  analysis.  Figure 
A-3.8  plots  N(.  versus  the  velocity  for  the  one-dimensional  analysis. 
Typically,  the  number  of  C-P  cycles  used  should  be  approximately 
one-half  the  value  of  N^.  for  the  maximum  velocity  anticipated  at 
the  measurement  location,  l.e.,  the  average  velocity  plus  3  tur¬ 
bulence  intensity  increments.  For  example,  if  the  average  velo¬ 
city  was  10  m/s  and  T.I.  =  505K  then  the  maximum  velocity 
anticipated  would  be  25  m/s.  From  Figure  A-3.8,  N^.  =  2^  at  25  m/s. 
Therefore  the  C-P  could  be  set  for  2^  cycles  and  the  uncertainty 
in  the  measurement  would  be  reduced  by  a  factor  of  4  from  that 
computed  previously. 

A-3.1.5.3  Data  Uncertainty  Due  to  Noise 

The  analysis  presented  above  has  assumed  that  the  S/N  was 
sufficiently  high  that  the  uncertainty  due  to  the  noise  in  the 
signal  was  negligible  compared  with  the  C-P  timing  uncertainty. 

Often  this  is  not  the  case.  The  S/N  degenerates  as  the  frequency 
increases.  Two  factors  determine  the  S/N.  First,  the  signal 
decreases  in  direct  proportion  to  the  frequency  increase. 

Second,  the  noise  amplitude  will  change  depending  upon  its 
origin.  If  the  major  noise  source  is  detector  shot  noise,  then 
the  noise  will  decrease  as  the  square  root  of  the  signal,  so  that 

S/N  a  (wd)-1/2 

On  the  other  hand  the  signal  may  decrease  to  the  point  that  ther¬ 
mal  noise  becomes  dominant.  As  the  Doppler  frequency  increases, 
the  amplifier's  bandwidth  has  to  increase  proportionally  and  the 


VN 


thermal  noise  increases  as  the  square  root  of  the  bandwidth.  In 
the  limit  of  the  thermal  noise  being  dominant 

S/N  a 

This  is  the  most  common  situation  in  aerodynamic  LDA  studies.  A 
study  of  the  effect  of  noise  on  the  LDA  frequency  determination 
has  been  made  by  Adrian  et  al.24  whiffen  et  al.25  have  investi¬ 
gated  a  model  proposed  by  Mayo  in  which  the  apparent  turbulence 
intensity  (En) ,  due  solely  to  the  noise,  can  be  expressed  in  terms 
of  the  S/N  for  a  C-P  using  N^  Doppler  cycles: 

En  =  - -  .  (A-35) 

4]f/SNR  N 

c 

Their  measured  results  agreed  within  +50%  of  the  prediction. 

Based  on  this  relation  a  minimum  S/N  of  4  is  required  to  ensure 
the  noise  uncertainty  contribution  to  the  velocity  is  less  than 
0.5%  when  2^  cycles  are  used  by  the  C-P.  Pfeifer^B  has  made 
similar  observations  in  a  study  in  which  he  also  examined  the 
influence  of  the  comparison  accuracy  for  a  double  counter  C-P 
system,  i.e.,  how  the  uncertainty  is  affected  by  the  accuracy 
setting  on  the  5:8  comparison.  His  results  show  that  the  uncer¬ 
tainty  is  independent  of  the  comparison  accuracy  and  he  measured 
uncertainties  comparable  to  those  of  Whiffen. 25 

A-3 . 1 . 5 . 4  Statistics  Uncertainties 

Beyond  the  uncertainties  in  a  single  measurement  there  are 
also  the  questions  of  how  many  measurements  are  required  to 
adequately  define  the  velocity  distribution  and  of  how  long  measure 
ments  must  be  taken  in  the  flow  to  fully  represent  all  possible 
flow  conditions  at  the  location.  The  question  of  the  number  of 
measurements  was  addressed  by  Cline  and  Bentley. 27  They  showed 
that  if  the  velocity  distribution  could  be  characterized  as  a 
one-dimensional  normal  distribution  with  mean  value  )  and  stan¬ 
dard  deviation  a  (turbulence  Intensity  TI  =  »i),  then  if  a  sample 
of  size  N  gave  a  mean  value  of  v  and  a  standard  deviation  of  s, 


the  variance  in  v  is 


af  =  a^/N  (A-36) 

V 

and  the  tolerance  in  v  is 

Tol(v)  =  Z  ^  =  Z  TI/Zn  .  (A-37) 

Thus  a  955K  assurance  that  |i  is  within  Tol(v)  around  v  is  achieved 
when  Z  =  1.96  (normal  distribution).  Figure  A-3.9  is  a  plot  of 
Equation  (A-37)  for  several  TI’s,  with  Z  =  1.96.  For  example,  to 
be  95%  sure  that  v  is  within  5%  of  m  at  a  measurement  location 
where  the  TI  is  20%  requires  more  than  50  random  measurements. 
Higher-order  moments  will  take  longer  to  stabilize.  The 
tolerance  in  the  turbulence  intensity  is  expressed  as 

Tol.  (TI)  =  Z  / +  (TI)^]  ,  (A-38) 

where  /3  is  the  kurtosls  of  the  distribution  (the  fourth  moment 

about  the  mean  divided  by  the  square  of  the  second  moment  about 

the  mean).  The  assumption  that  <0.3  has  been  made  to  obtain 

n 

Equation  (A-27).  For  a  normal  distribution  b  *  3.  If  the  distri¬ 
bution  is  not  normal  the  implication  of  the  value  of  Z  must  be 
re-evaluated.  Assuming  a  normal  distribution.  Equation  (A-38)  is 
plotted  in  Figure  A-3.10.  To  be  95%  assured  that  the  measured  TI  is 
within  5%  of  the  true  TI  requires  more  than  500  random  measure¬ 
ments,  an  increase  by  a  factor  of  10  over  the  requirements  for 
determining  v  to  the  same  accuracy  and  assuredness.  The  number 
of  measurements  will  Increase  further  if  higher  moments  are 
desired.  Also,  in  many  measurement  situations  the  distributions 
are  not  Gaussian  and  fi  can  be  substantially  larger  than  3.  As  a 
result,  the  requirement  for  the  number  of  measurements  increases 
still  more. 

A-3.1.5.5  Uncertainties  due  to  Flow  Characteristics 

The  preceding  discussion  has  assumed  that  the  velocity 


measurements  are  independent  (random).  Edwards^®  has  Investi¬ 
gated  the  effect  on  v  of  the  measurement  rate,  the  flow  proper¬ 
ties,  and  the  C-P  properties.  He  found  that  the  variance  in  the 

2 

estimate  of  the  measurement  of  the  mean  velocity,  a  ,  is  given  by 


2Rt  ^ 
c 


(A-39) 


where  is  the  flow  correlation  time,  T  is  the  time  over  which 

2 

measurements  were  made,  a  is  the  variance  in  individual 

m 

measurements  due  to  errors  (bias  or  uncertainty) ,  R  is  the  data 

2 

arrival  rate,  and  is  the  flow  variance.  The  key  to  accurately 

determining  the  average  velocity  is  to  take  measurements  for  many 
correlation  times,  T>>t^^  Merely  increasing  R  is  not  sufficient. 
The  combination  of  R  and  T  should  be  chosen  to  provide  sufficient 
data  elements  to  reduce  statistical  uncertainty  while  also 
reducing  the  uncertainty  in  the  mean,  due  to  flow  conditions, 
both  to  within  acceptable  bounds.  In  general,  Edwards  suggested 
that  the  arrival  rate  should  be  adjusted  to  ensure  that  T  >  200 
Then,  based  upon  the  computed  value  of  T(.  and  several  parameters 
characteristic  of  the  apparatus,  a  corrected  v  could  be  computed. 

The  key  point  that  should  be  made  is  the  interactive  nature 
of  the  data-taking  process.  The  amount  of  data  needed  to 
correctly  characterize  the  flow  at  the  measurement  location 
depends  upon  the  flow  being  measured.  Consequently  an  interac¬ 
tive  computer  with  on-line  data  reduction  is  required  to  deter¬ 
mine  whether  the  measurement  made  meets  the  accuracy  desired. 

A-3.1.5.6  LDA  Beam  Deflection  Effects 


Another  potential  source  of  error  in  the  dual  beam  LDA 
system  would  be  the  motion  of  the  beams  due  to  density  fluc¬ 
tuations  along  the  optical  path.  This  will  also  give  rise  to  the 
wd  contribution  in  Equations  (A-29a)  and  (A-29b) .  Two  other  poten 
tlal  error  modes  have  received  some  study.  The  "fringe"  spacing 
could  change  altering  the  frequency  to  velocity  conversion  fac- 
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tor,  K  [Equation  (A-16)],  and/or  the  "fringes"  could  be  in  notion 
contributing  an  additional  frequency  offset.  Bouis^^  monitored 
the  "fringe"  system  in  a  hot  flow  supersonic  wind  tunnel  (1700*K, 
720  m/s)  and  found  that  the  "fringe"  spacing  fluctuation  was  f 
O.IE  and  the  combined  sources  of  error  gave  an  uncertainty  of 
0.3X.  Both  of  these  errors  are  small. 


A-3 . 1 . 5 . 7  Seed  Inf luence 

The  size  requirements  on  the  seed  particles,  to  ensure 
faithful  tracking  of  the  velocity  field,  follow  directly  from  an 
analysis  of  the  relative  motlon^O  of  the  fluid  and  the  seed  par¬ 
ticle.  In  gas  flows,  where  Pp/Pf  >>  1  the  equation  of  motion 
reduces  to 


P 


(A~40) 


where  p  is  the  density,  U  is  the  velocity,  and  d  is  the  diameter. 
The  subscripts  are  p  for  particle  and  f  for  fluid.  Quite  often 
the  criterion  used  to  evaluate  a  particle's  suitability  is  based 
upon  the  fidelity  of  the  particle's  kinetic  energy  compared  to 
that  of  the  fluid.  This  has  the  form 


0 


U. 


1+ 


2irf 


(A-41) 


where  f  is  the  turbulent  fluid  fluctuation  frequency  and 


C  = 

P  d 
P  P 


^  .  Typical  response  curves  are  shown  in  Figure  A-3. 11. 


It  should  be  understood  that  this  test  is  only  on  the  kinetic 
energy.  There  are  more  stringent  tests  that  require  fidelity  in 
both  amplitude  and  phase  for  the  particle's  response  to  the  fluid 
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fluctuations.  These  tests  and  their  results  in  subsonic  and 
supersonic  flows  have  been  revlewed^^-'^^  and,  as  a  rule  of  thumb, 
Melling^l  suggests  that  subsonic  flows  be  seeded  with  particles 
<1  (la  in  diameter  and  supersonic  flows  be  seeded  with  particles 
0.3  to  0.5  urn  diameter.  A  major  difficulty  is  the  decrease  in 
scattered  intensity  as  the  diameter  decreases.  For  the  diameter 
mentioned,  the  scattered  intensity  is  approximately  proportional 
to  d^  and  it  may  be  necessary  to  trade  response  for  signal 
depending  upon  the  circumstance.  Devices  for  dispersing  and 
injecting  the  seed  are  reviewed  by  Helling. 

Seeding  the  flow  appropriately  and  the  ability  of  the  seed 
to  follow  the  flow  have  been  major  concerns  in  wall  boundary 
layer  studies.  Many  studies  have  measured  a  decrease  in  the 
Reynolds  shear  stress  very  close  to  the  wall  when  using  LDA 
instrumentation.  This  is  contrary  to  prediction  and  to  hot  wire 
measurements.  Demotakis^^  has  suggested  that  this  error  is  seed 
related.  He  suggests  that  the  seeds  do  not  follow  the  strong 
turbulent  structures  near  the  surface.  Instead  they  impact  the 
surface  and  stick  and,  as  a  result,  the  LDA  measurements  are 
biased,  showing  a  net  velocity  into  the  wall.  These  results  have 
been  called  into  question  by  Schairer^^  who  suggests  that 
orienting  the  LDA  fringe  system  at  +45*  to  the  wall  normal  and 
using  a  Bragg  cell  frequency  offset  might  compensate  for  LDA 
measurement  errors  taken  in  the  boundary  layer  very  near  the  sur¬ 
face.  His  measurements  had  too  much  scatter  to  be  conclusive  but 
the  trend  appeared  correct.  This  orientation  and  frequency  off¬ 
set  of  the  LDA  system  are  also  suggested  in  this  report.  The 
result  may  be  doubly  beneficial. 

A-3.1.5.8  Scattering  Considerations 

Another  Important  consideration  in  the  design  of  the  LDA 
system  is  the  placement  of  the  detector.  Figures  A-3.12a-c  are  plots 
of  the  scattering  intensity  as  a  function  of  azimuthal  angle  for 
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COLLECTION  ANGLE 


Figure  A 


3.13.  Scattering  factor  as  a  function  of  collection 

angle  for  Silicone  Oil  droplets  of  0.5,  1.0,  and 
1.5  |im  diameter  and  a  detector  collector  centered 
at  165<*  to  the  forward  LDA  optic  axis. 


the  components  perpendicular  and  parallel  to  the  scattering 
plane. The  radial  distance  scales  logarithmically  and  It  can 
be  seen  that  there  are  large  variations  In  amplitude  over  small 
azimuthal  angle  variations  near  the  backward  direction.  The 
plots  are  for  silicone  oil  droplets  0.5,  1.0,  1.5  )m  In  diameter, 
and  Input  radiation  at  514.5  nm.  If  the  seed  Is  monodlspersed , 
the  detector  can  be  positioned  to  collect  a  maximum  In  the  scat¬ 
tered  radiation.  If  there  Is  a  spread  In  the  seed  diameters,  the 
detector  position  must  be  a  compromise  choice.  There  Is  always  a 
local  meoclmum  In  the  backward  direction,  but  using  this  lobe 
results  In  the  worst  case  for  LDA  spot  localization.  Collecting 
the  backseat tered  radiation  at  an  angle  may  result  In  the  detec¬ 
tor  being  placed  over  a  local  minimum  In  the  scattering 
amplitude,  such  as  at  150*  for  the  0.5  pm  droplets,  or  165*  for 
the  1.0  urn  droplets.  Therefore,  scattering  plots  should  be  used 
to  help  select  the  position  for  the  collector.  The  ^Information 
can  be  used  to  advantage.  Often  seeders  will  provide  droplets 
with  a  narrow  distribution  of  diameters.  If  the  smaller  scat- 
terers  will  follow  the  flow  better  than  the  larger  scatterers, 
the  collection  angle  can  be  chosen  to  enhance  the  signals  from 
the  desired  size  scatterer,  e.g.,  choosing  a  collection  at  165* 
discriminates  against  droplets  1.0  ixm  diameter.  But  If  the  range 
of  diameters  Is  too  large,  this  type  of  discrimination  Is  not 
effective  (Fig.  A-3.12c). 

It  should  also  be  remembered  that  the  scattered  Intensity 
is  strongly  dependent  upon  the  size.  Therefore,  even  though 
there  may  be  a  minimum  in  the  local  scattering  amplitude,  used  to 
discriminate  against  a  particular  diameter,  the  actual  scattering 
amplitude  for  the  size  discriminated  against  may  not  be  less  than 
for  the  desired  size.  Figure  A-3.13  Is  a  plot  of  the  collected  scat 
taring  Intensity  for  different  sized  collection  angles  centered 
at  165*.  Although  this  angle  represents  a  minimum  for  the  1.0  pm 
droplets,  the  scattered  Intensity  for  the  two  diameters  Is  almost 
the  same.  As  the  collection  aperture  increases,  so  that  the 


collected  scattered  radiation  includes  that  from  regions  beyond 
the  minimum,  the  scattering  intensity  from  the  1.0  (im  diameter 
droplet  increases  much  more  rapidly  than  that  from  the  0.5  tim 
droplet  and  becomes  an  order  of  magnitude  larger  for  a  collection 
half-angle  of  10* . 

Another  parameter  of  Importance  is  the  signal  quality. 
Positioning  of  the  detector  and  selection  of  the  scatterer 
diameter  should  not  be  based  solely  on  the  scattering  amplitude 
but  should  also  take  into  account  the  Doppler  signal  amplitude, 
especially  by  looking  into  the  expected  S/N  at  the  detector. 

Shown  in  Figure  A-3.14  is  a  plot  of  the  Doppler  signal  visibility, 
the  pedestal  peak  rms  amplitude,  and  the  maximum  S/N  in  the  case 
of  shot  noise  limitation.  Figure  A-3.14  was  plotted  for  the  detec¬ 
tor  located  at  165*  azimuth.  It  is  evident  that  using  a  scat¬ 
terer  with  radius  0.35  )m  to  0.70  mb  results  in  a  significant 
decrease  in  the  visibility  and  will  show  a  corresponding  decrease 
in  the  S/N  compared  with  signals  from  scatterers  with  radii 
either  slightly  smaller  or  larger. 

The  factor  needed  to  convert  from  the  relative  scattering 
amplitude  scale  used  in  Figure  A-3.14  is  k“2  (jc  =  2n/X).  Therefore 
a  pedestal  amplitude  of  1  corresponds  to  a  scattering  cross- 
section  of  6®7  X  10~1®  m2  at  514.5  nm.  For  a  total  laser  power 
of  1  Watt  into  the  250  )m  diameter  spot,  the  average  number  of 
photons  collected  by  the  detector,  at  the  location  and  with  the 
dimensions  given,  would  be  3®53  x  10.^^  To  resolve  a  Doppler 
cycle  the  signal  should  only  be  integrated  for  1/6  of  a  cycle. 

At  70  MHz  this  computes  to  567  photons  per  integration  period. 

This  will  vary  depending  upon  whether  the  scatterer  is  at  the 
center  of  the  LDA  spot  or  near  the  edge  and  whether  the  maximum 
or  the  average  of  the  Doppler  cycle  is  being  considered. 
Nevertheless  the  computation  is  a  good  indicator  of  the  expected 
signal  strength  and  the  associated  Poisson  fluctuation.  A  more 
accurate  computation  will  depend  upon  the  particular  geometry 
used,  the  scatterer 's  properties,  and  the  laser  power  in  the 
actual  system. 


A-3.2  ALTERNATE  TECHNIQUES  FOR  FLOW  FIELD  VELOCITY  MEASUREMENTS 


The  tlme-of-f light  techniques  were  introduced  at  the 
beginning  of  this  section  with  reference  to  the  two-spot  and 
transit  velocimeters .  There  is  another  class  of  techniques  based 
upon  time-of-f light  type  measurements  but  it  offers  a  distinct 
advantage  to  the  aforementioned  methods,  principally  it  records  a 
two-dimensional  field  of  the  flow  at  one  time.  The  basic  tech¬ 
nique  is  to  record  a  "snapshot"  of  the  flow  at  two  instants  of 
time  and  to  measure  the  change  at  each  location  in  the  flow  that 
is  imaged.  There  are  three  interesting  variations  of  this  tech¬ 
nique  which  will  now  be  discussed. 

First,  and  perhaps  most  obvious,  is  to  use  sheet  light  scat¬ 
tered  from  seed  particles  in  the  flow.  The  light  is  pulsed  twice 
and  a  simple  frame  of  film  is  double  exposed. 3®  The  pulses  are 
timed  so  that  the  motion  of  the  scatterers  can  be  seen  but  it  is 
not  so  great  that  there  is  confusion  interpreting  where  particles 
moved  to.  If  the  flow  is  seeded  with  particles  having  the 
appropriate  characteristics  (size,  shape,  etc.)  this  technique 
combines  both  visualization  and  quantitative  measurement. 

Another  appealing  feature  of  this  technique,  besides  its  two- 
dimensional  mapping  (velocity  and  f ield-of-view) ,  is  that  it 
le^  's  itself  to  straightforward  machine  analysis  of  the  data.  A 
laser  scanner  can  be  adjusted  in  size  so  that  it  can  illuminate 
only  one  exposure  pair  at  a  time  and  the  two-dimensional  velocity 
can  be  read.  A  variation  which  we  suggest  is  to  provide  some 
depth  to  the  light  sheet  and  to  photograph  the  scattering  from 
two  angles.  By  triangulating  the  resulting  velocity  vectors  it 
should  be  possible  to  measure  velocity  in  three  dimensions. 

The  second  technique  is  to  use  front- illuminated  holographic 
recording  of  the  seed  particles  in  the  flow.  If  a  double-pulsed 
hologram  is  recorded,  the  full  volume  illuminated  can  be  examined 
in  the  reconstruction  and  the  motion  can  be  determined  at  each 
plane  by  adjusting  the  focus  of  the  data  inspection  optical 
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system. The  advantage  to  this  approach  Is  that  the  entire 
flovrfield  is  recorded  at  one  time.  The  system  required  to  examine 
the  data  and  interpret  the  measurements  is  complex  but  the  poten¬ 
tial  is  so  great  that  it  is  being  pursued. 

The  last  technique  that  is  discussed  here  combines  some  of 
the  features  of  the  previous  two  but  with  a  major  distinction,  no 
seed  particles  are  used.  Instead  variations  in  density  in  the 
flow  field  are  used  as  velocity  tracers.  These  velocity 
variations  produce  an  effect  called  "speckle"  in  the  reconstructed 
wavefront  image. 38  In  a  double-pulsed  recording  the  density 
fluctuation  will  move  slightly  between  exposures  and  the 
reconstructed  speckle  pattern  will  exhibit  interference  fringes 
that  determine  the  local  two-dimensional  velocity. 

A- 3  SUMMARY 

The  various  techniques  for  measuring  velocity  have  been  sum¬ 
marized  in  Table  A-3.1.  The  range  of  operating  conditions  that 
each  method  can  be  used  to  study  is  so  great  that  only  the  coar¬ 
sest  of  scales  can  be  used  to  differentiate  the  techniques.  The 
principle  areas  of  concern  are  the  number  of  velocity  components 
determinable  and  a  combination  of  maximum  velocity  and  fractional 
velocity  that  need  to  be  measurable.  For  example,  hypersonic 
velocities  almost  require  a  standard  spectroscopic  technique  to 
measure  the  Doppler  shift.  But  in  these  cases  small  fractional 
changes  (<1%)  are  almost  impossible  to  measure.  The  overlap 
between  the  capabilities  is  quite  large  and  choice  of  a  par¬ 
ticular  method  will  depend  upon  the  total  system  requirements. 

In  only  one  situation  is  the  choice  clear  at  this  time — 
simultaneous  three  component  measurements  are  the  province  of 
LDA.  One  closing  comment.  Several  techniques  (CARS/SRS)  do  pro¬ 
vide  other  information  beside  velocity,  such  as  density, 
pressure,  and  temperature.  These  parameters  are  of  vital  concern 
in  aerodynamic  studies  and,  as  a  result,  some  combination  of 
techniques  may  be  highly  desirable. 


SECTION  A-4 
CONCLUSIONS 


As  a  result  of  our  evaluation  of  the  present  state-of-the- 
art  in  technologies  required  for  a  real-time  inter feranemography 
system,  we  have  found  that  although  we  can  define  the  laser 
Doppler  anemometry  system  for  Inclusion  in  real-time  inter- 
feranemography ,  it  is  not  possible  to  design  and  build,  at  this 
time,  an  interferometric  system  suitable  for  real-time  use  for 
all  the  types  of  air  flow  to  be  studied. 


A- 4 . 1  I NTERPEROMETRY 

Summarizing  the  prospects  of  obtaining  perturbation  values 
of  density  for  air  flows  in  wind  tunnel  applications,  the 
following  points  are  made: 

(1)  To  obtain  truly  accurate  perturbation  densities, 
multiple  views  of  the  field  are  required.  Without  multiple 
viewing  capabilities,  the  flow  fields  studied  have  to  be  no  more 
complex  than  two-dimensional  axisymmetr ic .  Even  then,  the  per¬ 
turbations  can  not  be  local  to  any  specific  point  in  the  flow; 
they  too  must  be  two-dimensional. 

(2)  In  addition  to  obtaining  multiple  views,  continuous 
interferometric  surveillance  of  the  flow  is  required.  Meeting 
this  requirement  using  holographic  interferometry  would  be  dif¬ 
ficult  with  its  requirements  for  vibration  Isolation,  separated 
two  beam  optical  geometry,  multiple  hologram  recording  and  pro¬ 
cessing,  and  temporal  coherence  requirements. 

(3)  The  use  of  Point  Diffraction  Interferometry  appears 
promising  as  a  technique  that  might  provide  continuous  inter¬ 
ferometric  surveillance  of  air  streams.  However,  this  technique 
has  not  been  applied  to  wind  tunnel  studies,  and  both  theoretical 
and  experimental  development  of  specific  PDI  processes  has  to  be 


accomplished  before  the  worthiness  of  this  technique  can  be 
assessed.  Therefore,  as  a  recommendation,  PDI  techniques  should 
be  analyzed  and  applied  to  control  experiments  to  determine  the 
degree  of  measuring  ability  especially  with  respect  to  sensing 
perturbations  in  density. 

(4)  Heterodyne  interferometry  techniques  offer  promise  to 
obtain  true  real-time  evaluation  of  the  object  wavefront. 

However,  this  technique  should  be  analyzed  and  applied  to  control 
experiments  to  determine  the  degree  of  real-time  measuring 
ability  especially  with  respect  to  sensing  perturbations  in 
density. 

(5)  Because  of  the  relatively  low  Interferometric  sen¬ 
sitivities  of  the  hypersonic  facilities  of  AFWAL/FI,  inter¬ 
ferometry  should  not  be  used  as  a  means  of  measuring  densities  in 
these  wind  tunnels.  Furthermore,  applications  for  subsonic  stu¬ 
dies  in  the  TGF  wind  tunnel  and  applications  at  the  high  total 
pressure  levels  of  the  Mach  3  and  6  wind  tunnels  also  appear 
poor.  The  best  applications  appear  to  be  for  transonic  and 
superso.iic  operation  of  the  TGF  facility,  and  for  moderately  low 
total  pressure  levels  of  the  Mach  3  and  6  wind  tunnels. 

A-4.2  LASER  DOPPLER  ANEMOMETRY 

This  study  has  investigated  the  requirements  for  making 
three-dimensional  LDA  measurements  in  a  wind  tunnel.  Where 
optical  access  is  limited,  as  in  the  AFHAL/FI  transonic  facility, 
the  optical  configuration  providing  the  best  measurement  accuracy 
is  shown  in  Figure  A-3.5a.  The  angle  between  the  LDA  systems, 
is  set  to  the  meucimum  allowable  consistent  with  good  experimental 
design.  Furthermore,  the  two-dimensional  LDA  system  should  be 
aligned  to  take  measurements  along  directions  at  ±45*  to  the  free 
stream  velocity.  In  boundary  layer  studies  it  is  assumed  that 
the  boundary  is  parallel  to  the  free  stream.  If  this  is  not  the 
case,  a  compromise  angle  has  to  be  used  that  is  consistent  with 


the  needs  of  both  the  free  stream  and  the  boundary  layer,  par¬ 
ticularly  when  near  the  boundary.  This  last  decision  is  based 
upon  the  requirements  for  accurate  Reynolds  stress  measurements 
near  the  boundary  and  has  been  discussed  earlier. 

It  has  been  shown  that  the  data  collection  and  data 
interpretation  are  Interactive.  Consequently  it  is  essential 
that  a  computer,  preferably  a  high-speed  system  capable  of  quasi 
real-time  data  evaluation,  be  used  for  data  collection  and  analy¬ 
sis.  In  addition  to  the  velocity  record,  other  parameters  should 
be  recorded  simultaneously.  These  include  the  time  at  which  the 
velocity  datum  became  available  and  the  number  of  Doppler  cycles 
in  the  burst.  Due  to  the  limited  capability  of  the  counter¬ 
processors  for  measurements  in  supersonic  free  streams,  alter¬ 
native  forms  of  processing  the  LDA  signals  should  be 
investigated. 
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In  order  to  determine  the  fluid  velocities  along  the  coor¬ 
dinates  of  interest  (x,  y,  z)  the  measurement  coordinates  frame 
(x' ,  y' ,  z')  must  be  known  accurately.  Often,  to  simplify  the 
analysis,  the  2-D  LDA  is  aligned  so  that  (x* ,  y' )  lie  in  the 
(x,  y)  plane.  In  the  LDA  configuration  where  a  non-orthogonal 
1-D  system  is  used  to  measure  will  be  the  largest  uncer¬ 

tainty.  By  ensuring  that  is  not  required  to  transform  V^' , 

Vy'  to  Vx,  Vy  the  uncertainties  in  the  result  are  minimized. 

When  the  x' ,  y'  £uces  are  aligned  with  x,  y  (Figure  A- la)  the 
situation  is  trivial.  When  the  flow  conditions  require  x' ,  y'  to 
be  at  +45*  to  x,  y  the  situation  is  not  so  simple  for  the  3-beam 
system.  Most  LDA  systems  effect  the  +45*  measurement  by  rotating 
the  beams,  as  shown  in  Figure  A-lb,  by  45*  around  the  optic  axis. 
Unfortunately  this  mixes  Vj.  into  the  equations.  If  the  crossing 
angle  is  small,  the  contribution  of  V2  to  and  Vy  can  be 
neglected,  but  this  is  not  always  true.  The  incident  beams  can 
be  positioned  so  that  the  contribution  of  V^'  to  and  Vy  is 
zero . 


Consider  the  optical  configuration  shown  in  Figure  A-2. 

The  beams  will  be  assumed  to  form  an  isosceles  triangle  at  the 
plane  of  the  lens  so  that  only  one  component  need  be  considered. 
V}('  is  measured  by  the  Doppler  signal  generated  by  the  beams  at 
(0,  yj-D,  o)  and  (xj,  yj ,  o)  due  to  scattering  occurring  at  the 
measurement  spot  (o,  o,  zj).  The  equation  of  the  plane  through 
these  three  points  is 
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in  the  reference  coordinate  frame.  The  direction  cosines  of  the 
normal  to  this  plaune  are 
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(A. 2) 


where  is  a  proportionality  constant.  The  measurement  direc¬ 
tion  x'  will  be  perpendicular  to  (A. 2)  and  also  perpendicular 
to  the  vector  that  bisects  the  intersection  plane  of  the  incident 

*1  D 

beams.  This  bisector  extends  through  — ,  Vj  "  2'  ® 


(0,  0,  z^)  and  has  direction  cosines 


(A. 3) 


The  direction  cosines  of  the  normal  to  (A. 2}  and  to  (A. 3)  are 
given  by  the  determinant 
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Therefore,  the  contribution  from  V  to  V  ,  will  be  zero  if 

Z  X 
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The  most  common  3-beam  configuration  uses  two  identical 
beamsplitters,  the  first  one  to  provide  two  beams  vertically  (y) 
separated  by  2d,  and  the  second  to  split  the  bottom  beam  horizon¬ 
tally  (X)  by  2d.  (This  will  not  give  +45*  but  rather  +63.4*.  It 
common  use  is  a  result  of  the  cost  of  available  beamsplitters 
versus  the  cost  of  custom-made  beamsplitters.)  In  this  case 
X.  =  d  and  D  =  2d  and  according  to  (A. 5) 
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Another  advantage  of  this  configuration  for  boundary  layer 
studies  is  that  the  value  of  df ,  dictated  by  the  free  stream 
velocity,  can  be  smaller,  enhancing  the  measurement  capability  in 
the  direction  normal  to  the  surface  when  inside  the  boundary  layer. 
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